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RESUMEN 

Este trabajo es el resultado de varios estudios de puentes regulares de concreto. A traves de ellos 
se trata de presentar en torma sencilla los caminos a tomar para plantear el problema sfsmico de 
dichos puentes, partiendo de la modelacion en el piano, con la correspond iente caracterizacion de 
las propiedades dinamicas, hasta tratar las distintas metodologfas de solucion a la respuesta sfsmica 
del puente. 



ABSTRACT 

This work is the result of several studies of regular concrete bridges. Through those studies the 
aim is to present, in a simple way, the path to take in order to expose the seismic problems that 
affect these bridges, starting from a plane modeling, with their corresponding characterization of the 
dynamic properties, until trying different methodologies of solution regarding the response of a 
bridge. 



1 INTRODUCTION 



Sismos como el de San Fernando de 1971 y el de Loma Prieta de 1989 causaron perdidas 
importantes en los puentes, originando que los analisis del comportamiento del puente bajo posible 
amenaza sfsmica se incenti varan y como consecuencia de esto, numerosas metodologfas se han 
venido desarrollando para determinar el comportamiento sfsmico de puentes en diferentes areas. 
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Muchos de los aspectos involucrados en los nuevos codigos de puentes corresponden al resultado de 
trabajos experimentales y anah'ticos, desarrollados despues de estos sismos. Sin embargo, algunas 
de las normativas tienen su origen en estudios de edificaciones y es de destacar la diferencia que 
pueden, y de hecho presentan, los puentes en comparacion con las edificaciones en la respuesta ante 
movimientos sismicos del terreno. Posteriormente sismos como el de Northridge de 1994 y el de 
Kobe de 1995 han resaltado aun mas la importancia de profundizar en el estudio del 
comportamiento sismico de puentes. 

En este trabajo se presenta una forma sencilla de abarcar diferentes aspectos dentro del estudio 
del comportamiento sismico de puentes, tales como la modelacion a efectos dinamicos, la 
caracterizacion dinamica, es decir, la determination de la rigidez, la masa y el amortiguamiento, el 
calculo de los modos y periodos de vibration fundamental, y las diferentes metodologias de calculo 
de la respuesta dinamica. El estudio esta orientado a puentes regulares, esto es, aquellos que no 
presentan variaciones abruptas o inusuales de masa, rigidez o geometrfa a lo largo de la estructura, y 
no cuentan con grandes diferencias de estos parametros entre apoyos adyacentes, excepto en los 
estribos. 

Por ser un estudio enfocado a puentes con estas caracterfsticas, el objetivo es presentar los 
diferentes aspectos de la manera mas sencilla y aplicable con pequenas herramientas 
computacionales. Es asi como la modelacion es planteada en el piano, diferenciando el 
comportamiento en las dos direcciones del puente, longitudinal y transversal. La respuesta de los 
puentes se estudia bajo comportamiento lineal, sin desconocer la importancia que puede tener el 
desarrollar analisis no lineales en puentes que asf lo requieran (Elnashai, 1996). A su vez, el puente 
es considerado con sus componentes mas elementales, no incluyendo el efecto de la interaction 
suelo estructura, la cimentacion y el suelo de acercamiento de los estribos. Aspectos que dado el 
tipo de puentes pueden ser ignorados en un primer acercamiento a la solution dinamica, como una 
manera de simplificar el analisis. 

En la realization de este trabajo ha servido como gufa de desarrollo normativas como Seismic 
Design Guidelines for Highway Bridges del ATC-6 de 1986, Guide Specifications for Seismic 
Design of Highway Bridges del AASHTO de 1983, Codigo Colombiano de Construcciones Sismo 
Resistentes de 1884, Standard Specifications for Highway Bridges de la AASHTO de 1983 y 1994, 
y Analysis and Design of Reinforced Concrete Bridge Structures ACI 343 R de 1988. 



2 MODELACION SISMICA DE PUENTES 



El tipo de idealization que se utilice para determinar el modelo dinamico de un puente debe 
caracterizar lo mejor posible el comportamiento real del sismo, y a su vez debe permitir describir su 
respuesta ante el sismo. 

El calculo de la respuesta dinamica se efectua utilizando un modelo matematico, que describa 
cuantitativamente las propiedades de masa, rigidez y amortiguamiento del modelo dinamico. Si se 
admite la hipotesis de continuidad del material, y si se tienen en cuenta las diversas formas 
geometricas complejas que pueden tener la estructura, es obvio que definir las propiedades 
dinamicas mencionadas resulta complicado. Ademas, el calcular la respuesta dinamica de una 
estructura implica establecer dicha respuesta en cada uno de los puntos de la estructura, es decir en 
una cantidad infinita. 

Luego, para facilitar el calculo numerico de la respuesta se puede preestablecer un numero finito 
de puntos para el calculo de la respuesta. El numero de desplazamientos en las diferentes 
direcciones en cada uno de los puntos corresponden a los grados de libertad que intervienen en el 
modelo matematico. 

Cada uno de los grados de libertad describe simultaneamente propiedades de masa y de rigidez. 
Generalmente las propiedades de rigidez se describen por medio de grados de libertad localizados 
en los nodos donde se interconectan los elementos o donde existen juntas, no coincidiendo con 
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la estructura, pues se estan describiendo con el mismo grado de libertad propiedades totalmente 
diferentes. Como solucion a lo anterior se puede trasladar las propiedades de masa a los nodos 
donde se interconectan los elementos (Garcia, 1980). 

2.1 Ecuacion de Movimiento 

Para analizar el comportamiento sismico, se inicia con el estudio de la ecuacion de movimiento 
de un sistema muy sencillo, de un grado de libertad, para despues ver lo que ha de suceder al 
modelar la estructura con varios grados de libertad. 

Sistemas de un grado de libertad. Un sistema estructural simple se puede idealizar como una 
estructura de un grado de libertad, siendo un ejemplo el comportamiento longitudinal de un modulo 
de un puente, el cual se puede idealizar actuando independiente del resto de la estructura, ver figura 

En esta ideal izacion se asume el diafragma, punto donde se concentra la masa, como un cuerpo 
rfgido, mientras las columnas son flexibles a la deformacion lateral pero rigidas en la direccion 
vertical. 

Al existir un movimiento de la estructura de un grado de libertad debida a una excitacion 
dinamica, en este caso producida por el movimiento inducido por un sismo sobre la base de la 
estructura, el desplazamiento importante seria el horizontal, u, de la masa respecto a la cimentacion. 
En el momento que este desplazamiento sea diferente de cero, se presentan fuerzas de inercia, 
fuerzas de amortiguamiento y las fuerzas mecanicas correspondientes a las relaciones esfuerzo- 
deformacion del elemento estructural, en el caso la pila. La fuerza de inercia depende de la 
aceleracion absoluta de la masa, mientras que las fuerzas de amortiguamiento y las fuerzas 
mecanicas dependen respectivamente de la velocidad y el desplazamiento relativo al apoyo. Al 
plantear el equilibrio de estas fuerzas se establece la siguiente ecuacion: 

m u + cu + k u = -mu g {t) (i) 

Donde u es el desplazamiento relativo, ii la velocidad relativa, ii la aceleracion relativa, 
tn la masa del sistema, c el coeficiente de amortiguamiento, k la rigidez lateral de la estructura. 

Sistemas de varios grados de libertad. Al estudiar los puentes como un sistema de varios 
grados de libertad, su idealizacion, como se muestra en la figura 2, se ve mas claramente al analizar 
la estructura bajo vibracion transversal. La masa concentrada a nivel del tablero esta denotada por 
nij, m 2t ...m n > donde m ; es la masa correspondiente al nodo j> nodos determinados de acuerdo a la 
modelacion deseada del puente. Las propiedades de rigidez de la estructura se caracterizan por 
rigideces laterales k lt k 2 ,... k n , donde k } es la rigidez lateral del nodo j (Chopra, 1980). La excitacion 
en este caso es el movimiento inducido por el sismo en la base de la estructura. Suponiendo que el 
movimiento en el terreno solo producira componentes horizontales. 

La ecuacion de movimiento puede ser escrita en forma matricial como: 



[m][u\ + [c][u\ + [/c][w] = -M[y][^(r)] (2 ) 

Donde [u\ es el desplazamiento relativo, [m] la velocidad relativa, [u] la aceleracion relativa, 
[m] la masa del sistema, [c] el coeficiente de amortiguamiento, [/:] la rigidez lateral de la 
estructura, [y] el vector de valores unitarios que definen la direccion de aplicacion del sismo. 
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SISTEMA ESTRUCTURAL CONSIDERADO 
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IDEALIZAClCN ESTRUCTURAL 







u g (t) 



P0SICI6N DEFORMADA DEL SISTEMA 



m (U t +Ug(t)) 



ku- 



mu+cu+ku==-mug (t) 
DIAGRAMA DE CUERPO LIBRE 



Figura 1. Sistema de un grado de libertad 
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IDEALIZAClON ESTRUCTURAL 



posiciOn deformada del sistema 



mi(u,+u 9 (t)) m 2 (u 2 +u 9 (t)) 




mU+cu+ku=-mu g (t) 
DIAGRAMA DE CUERPO UBRE 



Figura 2. Sistema de varios grados de libertad 
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3 MODELACION EN EL PLANO 

Se debe resaltar, que dada la posible ocurrencia de la excitacion sfsmica en cualquier direccion 
respecto a los ejes del puente, es necesario considerar la vibracion y sus efectos en las diversas 
direcciones de este. Como mfnimo deben considerarse la vibracion longitudinal y la vibracion 
transversal, las cuales son claramente identificables en puentes rectos y las que se defininan como la 
direccion de la cuerda que une los estribos, como la longitudinal, y la direccion normal a esta, la 
transversal en el caso de puentes curvos. 

Este caso mfnimo ocurre en puentes con geometrfa y distribucion de masa y rigidez regulares, 
para los cuales se pueden esperar un desplazamiento de la vibracion en las dos direcciones 
ortogonales, es decir una independencia que permitan vibrar longitudinalmente sin que haya 
movimiento transversal. Y el caso contrario, ocurre en puentes irregulares, donde la vibracion 
induce acoplamiento entre las diferentes direcciones. 

3.1 Comportamiento longitudinal 

En puentes rectos y en menor medida en puentes curvos, la vibracion longitudinal de la 
estructura esta relacionada con la consideracion del tablero como un diafragma que interconecta los 
elementos estructurales en los apoyos. El comportamiento longitudinal y la distribucion de fuerzas 
dependen de la continuidad de la losa del tablero y/o de la ubicacion de las juntas determinadas en 
el diseno. El diafragma se encarga de hacer que los elementos por el conectados operen como un 
conjunto; pudiendose plantear para el caso dos sistemas ilustrados en la figura 3. El primero 
corresponde a un modelo donde existen juntas construed vas transversales. La junta se modela como 
una articulacion, por lo cual se permite igualdad en los desplazamientos horizontales y verticales y, 
rotaciones diferentes en cada uno de los extremos de los elementos que conectan al nodo; 
impidiendose la transmision de momento entre vanos adyacentes. Los grados de libertad del 
modelo en cada nodo son tres, dos translacionales y uno rotacional, excepto para los nodos donde 
hay articulaciones, donde ademas del translacional cada extremo de elemento cuenta con un grado 
de libertad rotacional. 

El segundo ejemplo se muestra a la derecha de la figura 3, un modelo de sistema de tablero 
continuo que cuenta con estribo no monolftico. Con columnas empotradas en sus extremos. 

Las propiedades dinamicas de cada uno de los dos sistemas en manera muy resumida son 
presentadas a continuacion: 

Propiedades de Rigidez. La idealizacion estructural para el comportamiento longitudinal debe 
tomar en cuenta la rigidez aplicable de los elementos considerados. El tablero aporta las 
caractensticas de diafragma, por otro lado, de acuerdo a su rigidez relativa respecto a las pi las, 
puede generar un comportamiento de estructura de cortante, empotrando las pilas, o en el otro 
extremo no generar ninguna restriccion a la rotacion de la cabeza de ellas, con lo que estas se 
comportan como elementos sin continuidad longitudinal de momentos. 

De acuerdo al tipo de restriccion que se asuma que la losa le hace a las columnas, debera 
considerarse para cada una de ellas la rigidez. Para el caso doblemente empotrado (72£//L ? ), el caso 
con apoyo libre de rotar arriba ( 3EI/L *) o algun intermedio, siendo E el modulo de elasticidad del 
material, / la inercia y L la longitud del elemento. La rigidez total del conjunto sera la suma de las 
rigideces de los elementos que tienen desplazamientos comun a causa del diafragma. 
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Figura 3. Ejemplos de puentes modelados en sentido longitudinal y transversal 
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En los casos donde se presentan juntas en la losa del tablero, los tramos que la contienen aportan 
rigideces correspondientes a su situacion de elementos empotrados en un extremo y articulados en 
el otro. 

De existir apoyo continuo en los estribos, su consideracion implica tomar en cuenta la rigidez 
estructural del estribo. Por otro lado una idealizacion en la que no se considera el aporte del estribo 
a la rigidez longitudinal puede estar del lado seguro para el diseno de las pilas. En este caso, sin 
embargo, no puede olvidarse el efecto en el estribo, quien en realidad si va a recibir la carga 
longitudinal de existir el apoyo. Los resultados presentados en trabajos como los de Goel (1996) 
muestran que la rigidez de los estribos pueden ser significativamente diferentes durante las distintas 
etapas del sismo, mostrando un notorio decrecimiento a medida que las deformaciones aumentan. 

El empuje del terreno para la idealizacion longitudinal simple, puede ser ignorada. Sin embargo 
no se desconoce conclusiones como las expuestas en trabajos como los de McCallen et al. (1994) 
donde predicen las caracterfsticas de la vibracion natural del sistema de puentes indicando que el 
sistema natural de vibraciones es sensible tanto a la rigidez como a la inercia de los terraplenes del 
suelo. 

Propiedad de la Masa. La masa asociada al movimiento longitudinal es el correspondiente a la 
carga muerta de la fraccion del tablero que ha sido identificada como el diafragma, sumando parte 
de la masa de la infraestructura correspondiente. 

No es comun considerar la masa asociada a la carga viva, esto solo se haria en el caso de grandes 
proyectos viales en los cuales sea el caso comun contar con un uso pleno y continuo de la capacidad 
del puente (Bernal, 1985). 

Propiedad de amortiguamiento. Las propiedades de masa y de rigidez pueden ser calculadas 
en funcion de la forma y tamafio de los elementos estructurales, siendo imposible calcular de esta 
misma forma el amortiguamiento. En general al no ser posible el cuantificar la energia local de los 
mecanismos de disipacion, la matriz de amortiguamiento se expresa en funcion de la relacion de 
amortiguamiento, obtenidas por experiencias en estructuras similares. 

Para el caso de los puentes analizados se tomo un amortiguamiento del 5%. Trabajos como los 
realizados por McCallen et al. (1994) demuestran que niveles altos de amortiguamiento deben ser 
utilizados en un modelo sencillo lineal elastico con la finalidad de encontrar precision en la 
respuesta sismica del puente. 

3.2 Comportamiento transversal 

En la idealizacion estructural de un puente sometido a una vibracion transversal, el tablero es 
modelado como una viga apoyada sobre elementos flexibles, modelados como resortes. En estos 
resortes ubicados en los apoyos, se representa la contribution de la subestructura al desplazamiento 
transversal. De esta forma, se logra involucrar la continuidad de la superestructura y distribuir la 
carga sismica transversal en todos los elementos estructurales del puente de acuerdo con su rigidez. 
Los ejemplos idealizados para la vibracion longitudinal, son tambien modelados para la vibracion 
transversal. 

En el sistema continuo el comportamiento transversal es presentado en la figura 3, a la derecha. 
Los resortes se ubican en los puntos de conexion de pilas y estribos si estos se estan considerado. A 
cada nodo se le asignan tres grados de libertad, dos translacionales y uno rotacional. En la 
idealizacion se puede eliminar los grados de libertad en la direction longitudinal, despreciando de 
esta forma las deformaciones axiales en el tablero. 

En la figura 4 se muestra un modelo de ejemplo la idealizacion de elementos de la subestructura. 
La pila corresponde a un sistema de columnas empotradas a la viga cabezal y el estribo a un muro. 

El modelo de la izquierda en la figura 3, en los puntos donde se encuentran ubicadas las juntas, 
se colocan articulaciones. Es de escogencia arbitraria el numero de nodos al modelar el sistema. Sin 
embargo es de comun uso el considerar ademas de los puntos de conexion de elementos y puntos de 
juntas el adicionar modos cada tercio de luz, como minimo. 
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Las propiedades dinamicas de cada uno de los dos sistemas de manera muy resumida son 
presentadas a continuacion: 

Propiedades de Rigidez. La ideal izacion estructural para la vibracion transversal debe estimar 
la rigidez de todos los elementos, por lo cual la rigidez de los resortes corresponde a la rigidez al 
desplazamiento en la direccion transversal de las pilas y de los estribos en el caso que se tengan en 
cuenta. 

Esta rigidez puede determinarse empleando cualquier programa general de analisis matricial , 
que permita evaluar porticos. En el caso de estribos tipo muro, la rigidez puede ser determinada en 
la forma presentada en la figura 4 (Garcia, 1980). 

El tablero contribuye con las caracterfsticas de rigidez de una viga en cada tramo de luz. De 
existir apoyo continuo en los estribos, su consideracion implica tener en cuenta rigidez estructural 
del estribo. En los casos donde la juntas estan presentes, se deben tomar en consideracion para su 
aporte en la rigidez transversal, la condition de elemento empotrado en un extremo y articulado en 
el extremo opuesto. 

La falla de los diafragmas debilitan la rigidez transversal y la estabilidad de las vigas principales 
influyendo en el normal funcionamiento del puente. Los resultados, de la investigation de Xizhu, 
Chen y Rong (1996), indican que las fuerzas intemas y las deformaciones del diafragma son 
controladas por la vibracion a flexion transversal de las vigas principales del puente, y que una parte 
parcial de diafragma puede alcanzar la plasticidad mientras las pilas estan aun en estado elastico 
cuando se tiene una intensidad de sismo alta. Ademas, resalta como los efectos de las pilas sobre las 
fuerzas del diafragma principalmente pueden ocurrir en la vibracion transversal. 

Propiedades de Masa. La masa asociada a la vibracion transversal es la correspondiente a la 
carga muerta del tablero idealizado como viga. Al igual que en el analisis del movimiento 
longitudinal, el aporte por carga viva solo ha de tomarse para proyectos que asf lo justifiquen. 

Propiedad de amortiguamiento. Se le ha dado el mismo tratamiento que en la direccion 
longitudinal. 



4 EJEMPLO DE MODELACION EN EL PLANO 

Se presentan a continuacion un ejemplo de un puente vehicular, donde en forma descriptiva e 
ilustrativa se modela su comportamiento ante efectos sismicos. 

Ejemplo. Puente de hormigon armado para una via vehicular, sobre el Rio Fonce (Maldonado, 
1995), compuesto de seis porticos de hormigon armado configurados longitudinalmente, tres al 
inicio y tres al final, conectados por tres vigas presforzadas de section en I. Los porticos 
transversales se encuentran interconectados por una viga, ver figura 5. 

El si sterna estructural corresponde a un puente configurado por porticos y vigas preforzadas 
simplemente apoyadas a los extremos voladizos de las vigas de los porticos. La superestructura 
conformada por las vigas de los porticos, las tres vigas presforzadas y la losa. La infraestructura, la 
constituye las columnas de los porticos, los estribos en los extremos y sus correspond ientes sistemas 
de cimentacion. 

La modelacion longitudinal, corresponde a un portico piano, donde la viga lo conforman las 
vigas de los porticos, el tablero y las vigas presforzadas ; y las columnas, los elementos que 
constituyen las columnas de los porticos, ver figura 6a. 

Los puntos de conexion del diafragma con los estribos fueron modelados de diferentes maneras, 
como apoyos fijo-fijo, fijo-articulado, articulado-articulado e incorporation de los estribos. 
Finalmente los resultados mostraron (Maldonado, 1995) que el puente representaba un 
comportamiento mas real asumiendo la contribution de los estribos y tomando los puntos de 
conexion como nodos articulados. 

El puente en la modelacion transversal, trabaja la superestructura como una viga apoyada en 
elementos flexibles. Estos elementos flexibles se ubican en los puntos de conexion de estribos y las 
columnas, ver figura 6b. 
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Los porticos transversales formados por una viga y las tres columnas es modelada como un 
portico piano, y representan a cada uno de los resortes del modelo transversal, ver figura 6c. La 
rigidez de los estribos, a su vez, tambien fue tenida en cuenta en esta direction. 




Figura 5. Puente carreteable sobre el no Fonce (Colombia) 
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5 MODOS DE VIBRACION Y PERIODOS FUNDAMENTALES 



En la tabla 1 sc prcsenta un resumen de los periodos dc vibracion del ejemplo anteriormente 
planteado y de algunos mas realizados (Maldonado, 1991, 1995), modelados en el piano para la 
vibracion longitudinal y para la direccion transversal. En algunos casos, a su vez, son comparados 
con modelos tridimensionales. 

Los periodos en las dos direcciones fueron calculados siguiendo programas como CAL-86, 
SEISAB, o SAP90. El tipo de puente 1 se trabajo con SAP90 bajo un modelo tridimensional y con 
el metodo aproximado de la fuerza estaticamente equivalente; los resultados fueron muy similares. 
El puente 2 fue analizado con el programa CAL-86. Los puentes 3 y 4 con el metodo aproximado y 
los puentes 5 y 6 se trabajaron con las metodologias de analisis dinamico mediante CAL-86 y 
SEISAB respectivamente, y un segundo resultado de periodos se calculo utilizando la formulacion 
de Rayleigh (ATC-6, 1986). 

La formulacion de Rayleigh considera los siguientes aspectos: 
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1. Formulacion del modelo de parametros generalizados de un grado de libertad del puente en 
consideracion. 

2. Obtenida la modelacion de la estructura, se aplica una carga distribuida unitaria por unidad de 
longitud, P () y la deflexion resultante v s (x) es determinada. 

3. Se calculan a , /3,y por las expresiones: 



a = J v s (x)dx 


( 3 ) 


P = j w{x)v s (x)dx 


( 4 ) 


Y = J w ( j :) v ,(^) 2 dx 


( 5 ) 



Donde, w(x) es el peso por unidad de longitud, y los factores a , /J,y en unidades de longitud 

al cuadrado, fuerza por longitud y fuerza por longitud al cuadrado, respectivamente; de esta manera 
se llega al valor del periodo T: 

T=2K \f^~ »> 

\ P 0 ga 

Siendo, g la gravedad. 

En los valores de la tabla 1 se utilizo el espectro de diseno elastico del Codigo Colombiano de 
Construcciones Sismoresistentes (1984). La ultima columna de dicha tabla equivale al tipo de 
metodologia utilizada para hallar los valores mencionados. En el siguiente numeral se comentaran 
dichas metodologfas. 



6 PROCEDIMIENTOS DE ANALISIS SISMICO DE PUENTES 



Estudiado el problema de idealizacion del sistema ante el efecto sismico, definido el sismo de 
entrada en la modelacion (acelerogramas o espectros elasticos de diseno), se vuelve a la 
incertidumbre de la metodologia apropiada a utilizar para evaluar el comportamiento dinamico de 
un puente ante la accion de un terremoto. 

Diferentes metodologfas se pueden emplear dependiendo de la complejidad de la estructura y del 
grado de exactitud de la respuesta esperada. Algunas de estas metodologfas enunciadas se describen 
a continuacion: 

Analisis Dinamico con Sismos Reales o Sinteticos. Definidas las matrices de masa, rigidez y 
amortiguamiento de la estructura y la aceleracion del terreno debida al sismo, la respuesta de 
desplazamiento de la estructura puede estar determinada por la solucion a la ecuacion del 
movimiento. 

Para un sistema de nodos, esta ecuacion contiene n ecuaciones diferenciales con n 
desplazamiento desconocidos u/t), j = 1.2. ..n . De aquf que al resolverse simultaneamente se tenga 
que contar con metodos computacionales. 

Sin embargo la solucion simultanea se puede omitir utilizando metodos modales o metodos de 
superposicion de modos, “ Analisis Modal”. Estos metodos generalmente son aplicables a analisis 
de repuesta dinamica en estructuras complejas en el rango elastico, en particular para analisis de 
fuerza y desplazamientos de sistemas de varios grados de libertad producidos por vibraciones de 
intensidad media, los cuales causan respuestas esencialmente lineales. 

Otros metodos de solucion a la ecuacion de movimiento de un sistema dinamico, ademas del 
analisis modal el cual es el mas utilizado en problemas lineales, son analisis en el campo complejo 
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de la frecuencia y, la resolucion directa de las ecuaciones del movimiento mediante tecnicas de 
integration paso a paso (Barbat et al. 1994). 



SECCION 

LONGITUDINAL 



SECCION 

TRANSVERSAL 



Periodo (seq.) 



Long. 



METODOLOGIA 



Trans. 



1 -PUENTE VEHICULAR 



10m 









1 11.23 o 8.75 30 ,8.75 o 11.25 c 

SISTEMA PORTICO— V)GA PRESFORZADA 



0.58 



2.19 



ANAUSIS 

DINAMICO 



0.58 



2.20 



APROXIMADO 
F.E.E 



2-INTERCAMBIADOR 



5r 



n 

20 



10.5 



2.5 2.5 



1.19 



0.32 



ANAUSIS 

DINAMICO 



PRIMER TRAMO CURVO. SEGUNDO TRAMO CONTINUO SOBRE PILA AP0RT1CADA 



3-INTERCAMBIAMR 



8.6 , 



0.495 



0.565 



APROXIMADO 
F.E.E 



30 



30 



30 



V1GA— CAJON. PILA AISLADA DE SECCION VARIABLE 



4— INTERCAMBIADOR 




0.265 



0.40 



APROXIMADO 
F.E.E 



VIGA-CAJON. PILA AISLADA DE SECCION VARIABLE 



5— PUENTE PEATONAL 



0.33 



0.16 



ANAUSIS 

DINAMICO 



m-i 

10 10 

VIGA-LOSA. PILA AP0RT1CADA 



1 

2 



0.38 



0.165 



APROXIMADO 
FEE 



6— PUENTE VEHICULAR 



0.59 



0.34 



Jr 



ANAUSIS 

DINAMICO 



25 



25 



25 



0.59 



0.36 



APROXIMADO 
F.E.E 



VIGA-LOSA. PILA APORTICADA 



Tabla L Periodos fundamentales para algunos puentes tipicos 
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Analisis Dinamico con Espectros de Diseno. En los casos donde no se puedan utilizar las 
simplificaciones propias de los metodos estaticos equivalentes o cuando no se cuentan con registros 
de sismos bien sean reales o artificiales, es necesario recurrir a metodos dinamicos utilizando 
espectros de diseno elastico. Este es el caso, cuando se presenta la posibilidad de ocurrencia de dos 
o mas modos de vibracion con peso relativo similar o cuando la geometna es tal que existe 
acoplamiento entre las direcciones principales del movimiento. En estas situaciones es imposible 
establecer un desplazamiento unico que defina el movimiento de la estructura, lo cual hace 
necesario el utilizar varios modos. 

El metodo modal provee resultados exactos cuando se utilizan historias de respuestas dinamicas 
durante un sismo y resultados maximos estimativos al adoptarse a espectros de respuesta. 

Analisis con Fuerzas Estaticas Equivalentes. La contribucion de los diferentes modos de 
vibracion en los desplazamientos y fuerzas laterales dependen de varios factores, incluyendo la 
torma del espectro de respuesta, y sus modos de vibracion natural. Sin embargo, generalmente la 
contribucion del primer modo, el fundamental, es mayor en comparacion con los otros modos. 
Razon por la cual una buena aproximacion a la respuesta es utilizar las respuestas del primer modo, 
siendo equivalente esto al “ Analisis con fuerzas estaticas equivalentes”. Este metodo puede ser 
utilizado para movimientos sfsmicos longitudinales y transversales, y desde el punto de vista de 
dinamica estructural, se reduce a un problema estatico despues de que se introducen las fuerzas 
inerciales; pudiendo ser usado para puentes tanto continuos como discontinues regulares, que 
responden predominantemente en el primer modo de vibracion, con condiciones de fronteras en 
estribos y con pi las intermedias. 

La determinacion de la fuerza estatica equivalente se puede desarrollar a partir del calculo del 
periodo por la formulacion de Rayleigh, anteriormente expuesto. Con el periodo asf obtenido se 
evalua el coeficiente sismico, Cs, en los espectros de respuesta, y con el se determina la fuerza 
sfsmica estatica equivalente P o(x). La cual representa el primer modo de vibracion: 

BC 

f oW=y w(*) V s (x) (7) 

Luego se aplica la carga P fx), como una una fuerza estatica, sobre la estructura y se determinan 
las fuerzas y desplazamientos en cada miembro. 

En la tabla 1, la ultima columna corresponde a la metodologfa utilizada en cada puente para el 
calculo de su respuesta sfsmica. Los analisis dinamicos con modelos tridimensionales se realizaron 
con SAP90 y SEISAB en los puentes 1 y 6, y los restantes analisis dinamicos se desarrollaron con 
utilizando modelos en el piano, es el caso del puente 2 y 5. Para los restantes casos se 
empleo el metodo aproximado. 



7 CONCLUSIONES 

A continuation se resaltan algunos puntos o ideas que marcan en forma significativa el 
comportamiento sismico en el piano horizontal de puentes regulares, y los cuales, fueron deducidos 
del estudio anteriormente expuesto: 

• El tipo de idealizacion que se utilice para determinar el modelo dinamico de un puente, debe 
permitir describir lo mejor posible $u respuesta ante el sismo. El independizar la idealizacion 
ante el comportamiento longitudinal y transversal en puentes regulares es una simplificacion 
factible y send 11a de aplicar, que muestra un buen acercamiento a la respuesta bajo modelacion 
espacial. Un ejemplo el cual ilustra esta conclusion se ve en la tabla 1 con los puentes 1 y 6 y 
mas claramente en los resultados a nivel de fuerzas expuestas en el trabajo de Maldonado 
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• La cantidad de nodos a considerar en los modelos longitudinal y transversal esta en funcion de la 
exactitud deseada en la respuesta, y a su vez, en las condiciones mismas que el puente por su 
geometrfa y rigidez le impongan. Es asi, como en el modelo longitudinal no necesariamente solo 
se deben consideran nodos los puntos de interconexion de los elementos o la existencia de 
juntas. Entre ‘estos deben colocarse mas nodos con sus correspondientes valores de propiedades 
dinamicas. 

• Los apoyos elastomericos, bastante comunes en puentes regulares, pueden ser modelados como 
resortes con rigideces equivalentes a sus caracterfsticas (Mena, 1995), en este trabajo todos los 
apoyos en general fueron modelados como nodos especiales con condiciones de desplazamiento, 
dado las restricciones impuestas. 

• La formulacion de Rayleigh para la determinacion del penodo fundamental en los puentes 
regulares deja ver una buena aproximacion en comparacion con los valores obtenidos por otros 
metodos. Esto se muestra en tabla 1, donde los valores de periodos para los puentes 1, 3 y 4 
fueron calculados por CAL-86 y comparados con los obtenidos mediante la formulacion de 
Rayleigh. 

• Para puentes regulares la metodologfa de la Fuerza Estatica Equivalente ofrece una primera 
aproximacion a la respuesta ante el movimiento sismico. Si se desea conocer la respuesta en los 
modos de vibracion diferentes al fundamental el metodo modal utilizando espectros de diseno 
elastico ofrece una buena altemativa. 

• Una combinacion de fuerzas sfsmicas ortogonales es usada para considerar la incertidumbre de 
la accion del sismo y la simultaneidad de ocurrencia de las fuerzas en las dos direcciones 
horizontales ortogonales. Las fuerzas y momentos sismicos elasticos resultantes en las dos 
direcciones seran los maximos de la combinacion de dos casos de carga: la primera combinacion 
es considerando el 100% de los valores absolutos obtenidos en fuerzas y momentos en la 
direccion longitudinal y el 30% en la direccion transversal y, el segundo caso de carga es 
tomando el 30% de los valores obtenidos en la direccion longitudinal y el 100% en la direccion 
transversal. Los maximos obtenidos de los dos casos de carga seran las fuerzas y momento de 
diseno. Esta combinacion de carga propuesta por el ATC-6 se asumio como acertada, al dar 
valores cercanos a los obtenidos bajo un analisis tridimensional (Maldonado, 1991). 
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RESUMEN 

Este trabajo describe un modelo elastoplastico acoplado con un modelo de dano mecanico 
e higrometrico que permite simular numericamente materiales compuestos. El modelo se 
basa en la teoria de mezclas de substancias basicas y permite evaluar el comportamiento 
interdependiente de los distintos materiales que conforman el compuesto. La simulacidn 
numerica de pavimentos flexibles es urfa de las aplicaciones numericas de este modelo. Se 
presenta una correccidn a la teoria de mezclas que permite tener en cuenta el arido de 
tamano pequeno. Este modelo se ha utilizado para simular numericamente un ensayo tipico 
para proyecto de mezclas en pavimentos flexibles. 

ABSTRACT 

The paper describes a hygro-damage-plastic constitutive model appropriate for numerical 
simulation of matrix composite material. The road pavement material is the particular 
application of this model. The theory is based on the mixture of the basic substances of the 
composite and allows evaluating the inter-dependence between the constitutive models of the 
different compounding. A simple modification of the classical mixing theory is presented for 
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the evaluation of the bulk composite materials with aggregate of particles, as road pavement. 
Some examples of application are given. 

1 INTRODUCCION 

En este trabajo se presenta un modelo constitutivo basado en la teoria de mezclas para la 
simulation del comportamiento no-lineal de materiales de matriz compuesta, especialmente 
orientado a materiales de pavimentos flexibles, compuestos de matriz bituminosa y agregado 
petreo, sometidos a cargas y degradation por humedad. Tambien se incluye en la teoria 
de mezclas, una formulation basica para tratar la relation de aspecto de los granos del 
agregado petreo. La ventaja de la formulation que se describe a continuation, reside en la 
potencialidad de poder tratar el material compuesto corno la combination de los distintos 
materiales simples que lo componen y luego incluir esta estructura de trabajo dentro de la 
tecnica de los elementos finitos. 

A continuation se describe la teorfa de mezclas para un material degradable y plastifi- 
cable, bajo la influencia de efectos de humedad. Luego se presenta una modification de esta 
teorfa para tener en cuenta la influencia del tamano y forma de las partfculas que confor- 
man la matriz y por ultimo se presentan distintos resultados que derivan de su aplicaciOn a 
mezclas bituminosas para pavimentos. 

2 TEORIA DE MEZCLAS DE SUBSTANCIAS 

La teorfa de mezclas clasica esta basada en la combinacion e interaction de las sustancias 
bdsicas que conforman el compuesto (Green and Naghdi, 1965) 

(Trusdell and Toupin, 1960) (Ortiz and Popov, 1982) (Oiler et al., 1994). Se considera que 
en cada punto material en analisis, participan todas las sustancias componentes a la vez y 
cada una de ellas con su propia ley constitutiva y en la proportion de volumen asignada. 
Esto permite combinar materiales con comport amientos diferenciados, donde cada uno de 
ellos presenta un comportamiento evolutivo gobernado por sus propias leyes. 

Una caracterfstica b&sica de la teorfa de mezclas es que la deformation total en un punto 
6ij es la misma para t.odos los componentes del compuesto, presentando la siguiente ecuaciOn 
de cierre o compatibilidad: 



Esta hipOtesis vale rigurosamente sOlo si se aplica la teorfa de mezclas para materiales 
con comportamiento en paralelo, ej.: matriz con fibras largas, hormi-gOn armado, etc. En el 
caso de matriz con fibras cortas o matriz compuesta con agregados petrcos etc. es necesario 
definir otra ecuaciOn de cierre (Oiler et al., 1995). Otra alternativa para tratar materiales 
de matriz compuesta, es realizar una correction en las propiedades de cada componente y 
mantener la ecuaciOn de cierre (ec.l). 

Los materiales compuestos que cumplen con la ec.l, satisfacen tambien la condition 
bdsica de la aditividad de la energfa libre de sus componentes (Trusdell and Toupin, 1960). 
Esto es: 



6{j — (e^)} — {tij) 2 — ••• ( e *j)n 



( 1 ) 
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n 




donde $ c es la energia libre correspondiente a cada una de las ”n” sustancias componentes 
de la mezcla, k c = dh . y - (a p + d) c son la fraction de volumen y las variables inter- 
nas del componente c esuno , donde a p conforma el grupo de variables internas pldsticas , 
^ = Jo [dmec + dhj dt la variable interna de dano producida por la evolucibn del proceso 
mecanico e higrometrico y h es la variable libre que contempla el contenido de humedad me- 
dido en peso. Todas estas variables corresponden a cada uno de los materiales que forman 
el compuesto. 

La definicibn de la tension para el compuesto S zj se obtiene a partir de la formulacibn 
del modelo hiperelastico garantizando la condicibn de disipacibn del segundo Principio de 
la Termodinamica (Malvern, 1969). Esto es: 



c o 0* 

Sii=m m 



V 



n 



- kcTri c 



dVc 

9Eij 



n n 

= Y, k cm°c (S tj ) c = Y k c< (C ijkl ( h , d)) c (E e kt ) c 

c= 1 C— 1 

y el tensor tangente del compuesto, resulta: 



(«&<). 

donde Cf jkl representa el tensor constitutive tangente para el componente c“ 



( 3 ) 



( 4 ) 



2.1 MODELO CONSTITUTIVO DE UN SIMPLE MATERIAL 
COMPONENTE 

Cada sustaneia basica incorporada dentro de la teorfa de mezclas puede seguir cualquier ley 
constitutiva. El caso particular aquf tratado sigue una ley constitutiva elastoplastica con 
dano ( degrad acibn de rigidez), producido por efectos mecanicos y/o higrometricos. Esto 
significa que la evolucion de la tensibn depende de la variable libre de deformacidn E de la 
humedad del componente h y del grupo de variables internas pldsticas a p y de dario d. 

El modelo constitutive que aquf se describe tiene capacidad para simular el compor- 
tamiento do materiales ineldsticos. En este caso particular, se tratard aquellos materiales 
elastoplasticos (fenbmeno de acumulacibn tie deformacibn permanente) con dano (fenomeno 
de degradacibn de rigidez) bajo efectos de carga y humedad. 

A continuacibn, se presents la ley evolutiva tensibn-deformacibn para cada componente. 
Con el objetivo de una simplificacibn en los indices, se suprimira el sub-mdice ”c”. No 
obstante, adviertase que cada variable se refiere a la propiedad de una de las sustancias que 
forman el compuesto. 
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2.2 MODELO ELASTOPLASTICO CON DANO 

El modelo constitutive) que aquf se describe brevemente, corresponde al denominado ”Mod- 
elo de Dano Pliistico” , publicado por (Lubliner et al., 1989), ampliado en Oiler et al. (1990) 
con una posterior mejora en el acoplamiento entre el dano y la plasticidad, realizado por 
(Luccioni et al., 1996). El modelo se basa en la hipbtesis de aditividad de la energfa libre 
entre la parte elastica ® e y la plastica Esto es: 

V = V e (Eij, Efj.h, d) + V? (a?) = ^0 [EtjCijuEl] + * p K, h) (5) 

Considerando la segunda ley de la termodinamica (desigualdad de Clasius-Duhem 
(Malvern, 1969)), juntamente con la definition de energfa libre expresada en la ec.5, resulta 
la siguiente definicibn de la tensibn: 



f)\D e 

S tj = m°^-=Q jkl (d, h ) 



\ 

Eki ~ Eki “ a ki (^~ hp) 

\ 7 



(6) 



donde ra° es la densidad del material simple, E kl , Eki , son los tensores de deforma- 

cion eldstico, total, pl&stico, e higrometrico respectivamente, d im < d < 1 la variable de 
dano mecanica e higrometrica, que evoluciona desde un umbral inicial menor que la unidad, 
hasta la unidad. El umbral inicial involucra los estados previos no mecanicos, como por 
ejemplo el efecto de disgregacion por humedad que se tratara posteriormente. Por ultimo, 
Cij k i (d, h) = Q (d, h) Cf jkl es el tensor- constitutive del material, que depende de la vari- 
able de dano y del contenido de humedad, S (d, h) una funcion que depende del dano y 
la humedad y S t j el tensor de tensiones para un punto del component e. Para completar 
la presentacion del modelo constitutive, se define el modelo constitutivo elastopl&stico y el 
modelo de dano. 



2.2.1 Funcion de plasticidad y potencial plastico 

Se define como funcion umbral limite de plasticidad y como funcibn potencial las formas 
clasicas (Tresca, Von Mises, Mohr Coulomb, Drucker Prager, etc. (Lubliner, 1990)), ho- 
mogeneas y de primer grado en las tensiones. Esto es: 



F (Sij,a p ,h) = HSiJ-KiaP'h) (7) 

G(Sij) = g (Sij) = cte . (8) 

donde f (Sij) y g(Sij) son las funciones de tension equivalente y potencial equivalente, 
definidas en el espacio danado, K (a p , h ) es el umbral de resist encia plastica que sigue 
una regia de evolucibn interna diferente para cada tipo de material (endurecimiento, ab- 
landamiento, etc.), a p = /q a p dt las variables internas del modelo plastico y h el contenido 
de humedad. 



Un mo del o constitutivo elastoplastico acoplado con dario mec&nico e higrornetrico 



23 



2.2.2 Flujo pldstico y ley de evolucibn de las variables internas 

El flujo pldstico, se define siguiendo la forma que establece la plasticidad clasica y las 
variables internas plasticas, resultan a partir de la regia de flujo: 



£' P 



= A 



dG 

dSa 



a p = A H (Sij, a p ) =A (h p ). 



dG 

v ’ ij dSa 



( 9 ) 

( 10 ) 



donde A es el denominado factor de consistencia plastica, G es la funcion de potencial 
plastico y H una funcidn escalar de argumentos tensoriales definida para cada variable 
interna del problema. 



2.2.3 Funcion de dario 

Se define como funcibn hrnite de daiio: 



G D (S lj ,d) = S(S ij )-f c (d) (ii) 

donde S (Sij) es la funcion umbral de daiio, definido en el espacio danado, f c (d) es el umbral 
de resistencia de dano que sigue una ley de evolucidn interna diferente para cada tipo de 
material (endurecimiento, ablandamiento, etc.) y d = Jq d dt la variable interna del modelo 
de daiio, que incluye fenornenos rnecanicos e higrometricos. 

2.2.4 Ley de evolucidn de la parte mecanica de la variable interna de dano 

La ley de evolucion para la variable de dano mecanica, sigue la forma (Oliver et al., 1990): 

. 0G D 

dmec P> (12) 

donde /t= Elites el denominado factor de consistencia de daiio, G° es la funcidn 

umbral de daiio y = Cf jkl E e k[ es la tensibn en el material virgen. 

2.2.5 Ley de evolucibn de la parte higrometrica de.la variable de dano 

La evolucibn del dano por efectos higrometricos, depende de una definicibn apropiada que 
relaciona el fenomeno de humedad, el tiempo y la degradacion. En este caso particular, se 
puede escribir: 



dh- ( fh (h, h, ho j ^ 



(13) 



La definicion de daiio mecAnico ec.12 y daiio higrornetrico ec.13, deben cumplir ambas 
con la siguiente condicidn forzada: 
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doncle // t es una funcidn que relaciona el contenido de humedad h con la degradacidn hi- 
grometrica 0 < d h < 1 y donde h 0 es la maxima humedad que puede absorber el medio 
continuo. La inclusion del fendmeno higrometrico puede realizarse teniendo en cuenta una 
variable independiente del daho mecanico (Oiler and Ofiate, 1996). 

2.2.6 Ley constitutiva tangente 

La derivada temporal de la tension (ec.6) , resulta: 



OSi 



Sij dEki Ekl 



dSij 

dh 



dSa , ds 



h +-ZT d 

dd 



ij j^P 
JT &kl 



dE kl 



(15) 



Sij— 



dCijki (d, h) OCj jki (d, h) • 

dd cl dh 



■E e 



'kl 



+C,jki {d, h) ■ ( Eki - Ekl ~ E ki ) 



(16) 



Eki 



Sustituyendo en la ec.15 el tensor constitutive por su expresion Cijki (d, h) — (d, h) 

y operando, resulta la variacion temporal de la tensidn, 



Sij— C'ijkl (Ekl C ijkl Ekl 

Aplicando la condicidn de consistencia plastica y de daho F— 0,G = 0 se obtiene: 

v T . \T .h 



Sij= (Cf jk i) T Ekl ~ (ctjki) Eki 



(17) 



(18) 



Los tensores constitutivos tangente (Cf jkl ) r del probleraa mecanico ( C^ kl ) T correspondiente 
al problema de daho higrometrico resiulta: 



/ore dG dF fie 

° i jrsdSrs dSyrrn ° mnkl 



(cfjkl ) —Cijki ^ 2]<(h\ dG i A F C e -S 

d lt hu dStu ^ OSmn ^mriTS d 
r e 9G ( dF r° + 

j ^7* ^ ijrs dSrs w/nkl dh \ kl J J 

[E'ijklA — Cijki QK / u \ OG ■ dF f<e dG 

V ' ±^id^" li nud^^V^°rnnr S dSrs 



(19) 



(20) 



donde Cf ]kl es el tensor constitutivo tangente, que contiene la information de la degradation 
del material (Oiler, 1991) (Luccioni, 1993). 

1 donde: 

Cf jkl = SC* ** + £ {Tfj • S k i) 

Cijki = 9C&* + * {Tfj • Ski) - i {T t h j - Ski) 

rpD _ d* ( OG dS \ f>Q 
- dd y~ 9 S OSr» J rsi 3 

rrh _ da ( r h \~ l 
T U = dh [ Q ij) 
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3 MODIFICACION DE LA TEORIA DE MEZCLAS PARA MATRIZ CON 
AGREGADOS DE CORTA LONGITUD 

Los materiales de matriz cornpuesta, no cumplen exactaniente con la condicidn de compat- 
ibilidad expresada en la ec.l , debido a quo la transference de tensiones entre las partfculas 
y el res to de la matriz no es perfecta, siguiendo por lo tanto una distribucidn como la 
que se muestra en la Figura 1. Esta situacidn se produce por que no existe una longitud 
adecuada del refuerzo que permit e la transmision de los esfuerzos. Para ello es posible el 
plant eo de otra ecuacion de cierre de deformaciones (Oiler et al., 1995), o bien modificar las 
caracteristicas mecanicas de los componentes (Jayatilaka, 1979). 




Figuia 1: Mecanismo de transference de tensiones entre la matriz y un cornponente 

En todos los casos de matriz cornpuesta, existe una longitud de transference minima, 
que garantiza la perfecta compatibilidad entre matriz y el arido, a esta longitud se de- 
nomina longitud crftica / c . Cualquier arido cuya dimension fundamental sea inferior a 
esta inagnitud, no participa plenamente en los mecanismOs de transference de esfuerzo 
(Jayatilaka, 1979). Asf, la longitud crftica se obtiene a partir de la condicion de equilibrio 
entre el agregado y la matriz: 



<^max 

=f 7^r r (2i) 

donde se observa que la longitud depende de la maxima relation de resistencia axial (5 niax ) 
respecto de la tangencial (r max ), multiplicada por el radio medio del agregado (r). 

Estableciendo la distribution de tensiones en el cornponente arido (Figura 1), es posible 
obtener la siguiente ecuacidn diferencial que expresa el equilibrio del mismo segun el eje 
axial de mayor longitud (Jayatilaka, 1979). 
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(22) 



donde P es la fuerza maxima de interaction entre el agregado y la matriz, H una constante 
que depende del ordenamiento de los agregados, C c el modulo de Young del cornponente y 
A c la section transversal media del cornponente. La solution de la ecuacion diferencial 22 
expresa la fuerza P como: 



P = C\ sin (/3x) + C 2 cosh (0x) + C c A c Eu (23) 

donde C\ y C 2 son constantes, (3 es un coeficiente que depende de C c ,H y A s “ , Gin es 
el modulo elastico transversal del compuesto y r ' la distancia media entre las partfculas 
de agregado (ver Figura 2). Las dos constantes se obtienen de las condiciones de borde 
p = 0 para x = 0 y x = l. Una vez obtenidas las constantes de integration, se escribe la 
correspondiente tensidn como: 




Figura 2: Relaci6n de aspecto que se considera en el drido 



Si — C c Eu 



cosh 



1 - 



(9-fr) 



cosh 



VO < x < 



(24) 



Para incorporar esta contribution del cornponente agregado como un valor medio en de 
la teoria de mezclas, se establece su promedio como: 



5 1 



= y / Sdx = 

l JO 



C c Eii 



(25) 



2 a — I H = / £Lu±. 2 tt — 

I- Y ° cAc \l Cc ' A c>n('7r) 
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siendo C c e 1 modulo de Young promedio del arido ”c” 3 . 

La definicidn anterior para el modulo de Young, permite extender la energia libre 
mostrada en la ec.5 , a una expresion analoga, donde para cualquier componente agregado, 
se define como: 



con: 



* = ** (EivEfj, h,d) + W(cP) = [EtjC ijkl E e kl ] + ¥*(o»\ft) (26) 

\ tanh (f ) 



C ijki (d, h) — Cjjki (d, h) 



& 

2 



(27) 



A partir del segundo principio de la termodinamica, se deduce entonces la tensidn de la 
misma forma que se ha hecho con la ec.6: 



o _ 0^ e 

u _ m dE i:j ~ Cljkl 



tanh 



1 - 



(f) 



SL 

2 



(E kl - E p kl - 



(28) 



y la ecuacion constitutiva tangente del componente agregado, resulta como en la ec.18: 




De esta manera, se introduce la correccidn por transferencia de tensiones entre matriz 
y fibra, dentro del tensor constitutive del componente de agregado. 

La implement acidn del modelo constitutivo propuesto en un programa de elementos 
finitos no lineal sigue los pasos estandar para estos casos (Zienkiewicz and Taylor, 1994). 



4 PARAMETRIZACI6n DE LOS MATERIALES 

Con la finaiidad de estudiar el comport amiento de una mezcla asfaltica mediante la teoria de 
mezclas de sustancias , se presenta a continuacion, y a modo de ejemplo, la parametrizacidn 
de las propiedades de los materiales componentes de dicha mezcla asfeltica. 

4.1 CORRECCION DE LAS PROPIEDADES MECANICAS POR RELA- 
ci6n DE ASPECTO 

A continuacion se describe el procedimiento utilizado para la parametrizacidn de las tres 
mezclas asfdlticas que componen el firme en funcion de sus curvas granulometricas. Se 
denomina C-l al compuesto formado por betun mas la curva granulometrica utilizada para la 
capa de base de un firme, C-2 al compuesto utilizado para constituir la capa intermedia y C-3 
al que incorpora la capa de rodadura. Estos tres compuestos utilizan curvas granulometricas 
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con dos partes que se distinguen segun su comportamiento (ver Tabla 1): La primera parte 
de la curva, correspondiente a los diametros comprendidos entre 25mm. y 5 mm. que 
puede verse afectada por el daho higrometrico, perdiendo capacidad mecanica en forma de 
degradaci6n de la rigidez (C ljk i = S(d,h)C§ w J, donde d depende, entre otros fenomenos, 
de la evolucidn de d h , que es la variable de degradation por action de la humedad 14. La 
segunda parte de la curva , correspondiente a los diametros comprendidos entre 2,5mm. y 
el tamiz ciego que es insensible a la accibn de la humedad. 

Por razones de simplicidad en el calculo, se consideran sblo algunos de los diametros 
que participan en cada curva granulometrica como se muestra en la Tabla 1. De acuerdo 
a la forma, las propiedades mecanicas de los granos y su participacibn en el compuesto se 
deducen las propiedades corregidas de los aridos. 

En Tabla 2, se presenta la granulometria simplificada y el correspondiente mbdulo de 
Young corregido en funcibn de la relacibn de aspecto de los aridos y de la caracterfstica de la 
matriz (betun). Asf, el mbdulo de Young del agregado en estado natural es C a = 300000 ^ 
y corresponde a un material simple, homogeneo y continuo. En el caso que este material 
participe en forma de grano, pierde parte de esta propiedad por condiciones de adherencia 
y transmision de tension entre matriz y arido, llegando a C£ arr = 1559, 131 l^r, en tanto 
el mbdulo de Poisson se mantiene const ante en v = 0, 15. 

Tabla 1: Diametros de las curvas granulometricas 



djrimi.] 


% en peso que pasa 




C-l 


C-2 


C-3 


25 


90 


100 


100 


20 


78 


94 


100 


12,5 


56 


71 


82 


10 


52 


67 


65 


5 


36 


50 


41 


2,5 


25 


36 


14 


0,63 


13 


19 


8 


0,32 


8 


12 


17 


0,16 


6 


8 


6 


0,08 


4 


6 


5 



Tabla 2: Correccibn en funcibn de la relacibn de aspecto del arido 



Correcion por fibra corta Cb = 80, 8529 


L 


/ 

r 


r 


A 


c mat 


Vmat 


Cmat 


0 


t’o 


C C 


2 


1,05 


1.00 


0,7854 


80,375 


0.4 


28.705 


0,12525 


0.005197 


1559.1311 


1.25 


0,65625 


0.625 


0,30679 


80.375 


0.4 


28.705 


0,20014 


0.005197 


1559.1311 


0.5 


0.2625 


0.25 


0. 04908 


80.375 


0.4 


28.705 


0,50102 


0.005197 


1559.1311 


0,032 


0.0168 


0.016 


0.000201 


80.375 


0.4 


28.705 


7.82852 


0.005197 


1559.1311 


0,008 


0,0042 


0.004 


1 .2566e-5 


80.375 


0.4 


28.705 


31.3140 


0.005197 


1559.1311 


0,0001 


0.0000525 


0,00005 


1 ,9635e>9 


80.375 


0.4 


28,705 


2505.12 


0.005197 


1559.1311 
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5 CARACTERISTICAS MECANICAS DEL COMPUESTO Y DE CADA 
COMPONENTE 

El modulo elastico del betun a 60 °C, resulta de una parametrizacion con el compuesto, en 
este caso particular con el correspondiente ensayo Marshall. Esto es, por ejemplo para la 
mezcla C- 2 : 



C Marshall — ~r A ' a C^° rT 



15 °W = °’ 04C6 + °’ ^ = * 0 , 375^ (30) 

Se puede comprobar del <ibaco de Van der Poel (Yoder and Witczak, 1975), que para 
cargas de corta duration, se obtiene a partir del modulo elastico de la mezcla, el mismo 
resultado que en la ec.30 para el modulo elastico del betun. 

La tension umbral de inicio del comportamiento inelastico, resulta tambien a partir de 
una parametrizacidn con el ensayo Marshall. Como ejemplo, para la mezcla C- 2 , se acepta 
que cl inicio del comportamiento no-lineal de la mezcla Sh arshall = 20 ^, se debe a la 
plasticidad en el betun (vease la ligera perdida de linealidad que exhibe la curva carga 
desplazamiento de la mezcla C-2 en los primeros pasos del ensayo -figuras correspondientes 
a los ensayos Marshall). Luego de un cierto punto del ensayo, se observa un fuerte quiebre 
e.i la curva, instante en el que se inicia el comportamiento ineldstico del 4rido. Esto ocurre 
cu&ndo el betun ya exhibe una deforrnacidn excesiva, permitiendo el contacto de los aridos 
entre sf. Por esta razon, y apoyado en la condicidn de compatibilidad del compuesto, ec.l, 
se puede obtener el umbral de plastificacidn del betun para el instante en que se inicia el 
comportamiento no-lineal de la mezcla. Esto es: 

^ Marshall ~ ^ ^ Marshall — ^ 



E e = Y 500 = 0-0133 (deformation al lhnite elastico) (31) 

= E e Cb = 0.0133 • 80,375 = 1, 071— (tension umbral a 60°) (32) 

Queda entonces por determinar el umbral de dano del drido. Dado que esto ocurre 
avanzado el comportamiento plcistico del betun, luego del primer umbral de comportamiento 
inelastico de la mezcla, s 6 lo se debe exigir el cumplimiento de la siguiente desigualdad: 

S y a > E e Ca° rr = 0, 0133 • 1559, 13 = 20 , 7-% ( 33 ) 

cm 1 v ' 

correlacionando con el ensayo Marshall, resulta el siguiente valor, que cumple con la condi- 
tion anterior: 



Sl^ 30 



kg_ 
cm 2 



(34) 



( onfiimandose que el irido inicia su degradaci 6 n, ya avanzada la plasticidad en el betun. 
De forma analoga se procede con los compuestos C-l y C-3. 
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5 CARACTER1STICAS MECANICAS DEL COMPUESTO Y DE CADA 
COMPONENTE 

El mddulo elastico del betun a 60 °C, resulta de una parametrizacidn con el compuesto, en 
este caso particular con el correspondiente ensayo Marshall. Esto es, por ejemplo para la 
mezcla C-2: 

C Marshall = k b C b + k a C C ^ r 



1500 ^ = °> 04Cb + °’ mC ™ r = 80, 375^1 (30) 

Se puede comprobar del dbaco de Van der Poel (Yoder and Witczak, 1975), que para 
cargas de corta duracidn, se obtiene a partir del modulo elastico de la mezcla, el mismo 
resultado que en la ec.30 para el modulo elastico del betun. 

La tensidn umbral de inicio del comportamiento ineldstico, resulta tambien a partir de 
una parametrizacidn con el ensayo Marshall. Como ejemplo, para la mezcla C-2, se acepta 
que el inicio del comportamiento no-lineal de la mezcla S\ !arshall = 20^, se debe a la 
plasticidad en el betun (vease la ligera perdida de linealidad que exhibe la curva carga 
desplazamiento de la mezcla C-2 en los primeros pasos del ensayo -figuras correspondientes 
a los ensayos Marshall). Luego de un cierto punto del ensayo, se observa un fuerte quiebre 
e.i la curva, instante en el que se inicia el comportamiento inelastico del arido. Esto ocurre 
cudndo el betun ya exhibe una deformacibn excesiva, permitiendo el contacto de los aridos 
entre sf. Por esta razdn, y apoyado en la condicidn de compatibilidad del compuesto, ec.l, 
so puede obtener el umbral de plastificacion del betiin para el instante en que se inicia el 
comportamiento no-lineal de la mezcla. Esto es: 

qy 

° Marshall *~ / Marshall * 



E ~ 1500 = 0-0133 (deformacibn al h'mite eldstico) (31) 

= E e Cb = 0.0133 • 80, 375 = 1, 071^^ (tensidn umbral a 60°) (32) 

Queda entonces por determinar el umbral de dario del arido. Dado que esto ocurre 
avanzado el comportamiento pldstico del betdn, luego del primer umbral de comportamiento 
inelastico de la mezcla, sdlo se debe exigir el cumplimiento’de la siguiente desigualdad: 



s y > gegcorr = Q> ^33 . 155Qj ^ = ^ ? ±9 

cm 2 



(33) 



correlacionando con el ensayo Marshall, resulta el siguiente valor, que cumple con la condi- 
ci6n anterior: 



Sn - 30 



kg 



cm i 



(34) 



confirnidndose que el arido inicia su degradacibn, ya avanzada la plasticidad en el betdn. 
De forma analoga se procede con los compuestos C-l y C-3. 
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6 SIMULACI6N numerica 

Para la simulacibn numerica del comportamiento de la mezcla asf&ltica, se ha reproducido 
numericamente el ensayo Marshall que permite determinar la resistencia a la deformation 
plastica de mezclas bituminosas” . Este ensayo puede utilizarse tanto para el proyecto de 
mezclas en laboratorio como para el posterior control en obra de la misma. 

6.1 DESCRIPCION DEL ENSAYO MARSHALL Y SU SIMULAClON NUME- 
RICA 

El ensayo Marshall consiste en determinar la resistencia a la aplicacibn de una carga en una 
probeta cilfndrica de 4 pulgadas de di&metro (101,6 mm) y 2,5 pulgadas de altura (63,5 
mm) a una temperatura aproximada de 60 °C . El ensayo permite determinar dos valores: 

• Estabilidad : carga necesaria para producir el fallo de la probeta, 

• Deformacibn : disminucibn del diametro expresado en mm. que experimenta una 
probeta entre el comienzo de la carga y el instante de rotura. 

El ensayo permite graficar curvas Cargos - Desplazamientos con las cuales es posible 
determinar los valores de estabilidad y deformacibn definidos anteriormente (Norma NLT- 
159/86). 

Para la simulacibn numerica del ensayo se ha discretizado un circulo de 101,6mm. de 
diametro y dos sectores de circulo constituidos en acero, que simulan las mordazas de la 
prensa utilizada en el ensayo. La malla de elementos finitos utilizada posee 431 nodos y 776 
elementos finitos triangulares lineales (Figura 3). 



6.2 MATER1ALES UTILIZADOS EN LA SIMULAClON 

En la simulacibn numerica se ha utilizado la teoria de mezclas de sustancias basicas. La 
capa intermedia C-2 y la base C-l del firme, se han modelado como materiales compuestos 
con tres componentes basicos cada uno de ellos: 



6.2.1 Cornpuesto C-l 

Material C-l. a: corresponde al arido comprendido entre 20 y 5 mm. de la curva granu- 
lometrica C-l. Participa en un 58,62% y presenta un dano inicial en funcibn de la humedad 
y es inec&nicamente degradable. 

Material C-l.b: corresponde al drido comprendido entre 0,32 mm. y el tamiz ciego de 
la curva granulometrica C-l. Participa en un 32,97%, siendo degradable por problemas 
mecdnicos, pero insensible a la influencia de la humedad. 

Material C-l.c: corresponde al betun. Este material participa en un 8,41% del volumen de 
la mezcla y presenta un comportamiento elastoplastico. 



Ln modelo constitutive elastoplastico acoplado con daiio mecclnico e higrom£trico 
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Figura 3: Malla de elementos finitos del ensayo Marshall para mezclas asfdlticas. 



6.2.2 Compuesto C-2 

Material C-2. a: corresponde al arido coniprendido entre 20 y 5 mm. de la curva granu- 
lometrica C-2. Participa en un 45,02% y present a un dano inicial en funcion de la humedad 
y es mec&nicamente degradable. 

Material C-2.b: corresponde al arido comprendido entre 0,32 mm. y el tamiz ciego de 
la curva granulornetrica C-2. Participa en un 45,02%, siendo degradable por problemas 
mecanicos, pero insensible a la influencia de la humedad. 

Material C-2.c: corresponde al betun. Este material participa en un 9,96% del volumen de 
la mezcla y presenta un comportamiento elastoplastico. 

6.3 RESULTADOS OBTENIDOS 

El problema se resuelve en tension plana, permitiendo el flujo del material en el sentido axial. 
El casquete exterior del cilindro es de acero y el interior esta conformado por cada una de 
las tres mezclas. Entre el acero y la mezcla bituminosa, se considera deslizamiento, por 
friction, gracias al control de la tensibn tangencial entre ambos materiales. Los resultados 
obtenidos en las simulaciones numericas se presentan en la Figura 4, Figura 5, Figura 6, 
Figura 7 y Figura 8. En todas ellas se observa una disminucion del mbdulo de carga K = | 
a medida que aumenta el dano inicial. Sin embargo hasta valores del 40% de dano inicial 
en el material 1 no se observa una gran perdida de resistencia, situacibn que se manifiesta 
luego de superar los 2 mm. aproximadamente de desplazamiento. 
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Se observa, en la Figura 4 y Figura 5 para la mezcla C-2 y en la Figura 6 y Figura / para 
la C-l, una p6rdida de rigidez y resistencia en funcidn de la carga y del dano producido por 
efectos higrom6tricos, tendiendo en ambos casos asint6ticamente a una capacidad soporte 
del 42% de la maxima en el caso del material virgen. Es necesario hacer notar, que la 
respuesta experimental coincide exactamente con las respuestas ’ sin dano” mostradas en 
las Figuras 4 a 7. 



Ensayo Marshall C-2 




Figura 4: Curvas de respuesta del ensayo Marshall para la mezcla C-2 0 < dh < 0, 5 

En la Figura 8 se observa la evolucidn de la capacidad portante en funci6n del dano 
inicial higromfetrico, medidas todas para un desplazamiento de 1,86mm. La elecci6n de 
este valor corresponde al desplazamiento donde se produce la maxima carga de los ensayos 
Marshall C-l, C-2. 

Existe una correlaci6n entre el dano inicial, el contenido de humedad y el tiempo en 
que se mantiene esta situacidn. Para ello es necesario realizar estudios de laboratorio que 
establezcan la relacidn d= ( f h (h,h,ho)} donde h es el porcentaje en peso del contenido 
de humedad de un punto del drido. 

En la 9-1 se observa la malla de elementos finitos y los materiales empleados en la simu- 
laci6n numdrica. En la 9-2, 9-3 y 9-4 se observa para la mezcla C-2 un mapa de la evoluci6n 
del dano en la mezcla para tres estados distintos del proceso de carga. Particularmente la 
-2 muestra el dano mecdnico en la mezcla C-2 partiendo de un dano higromdtrico inicial 
nulo. La escala que se observa muestra el dano mecdnico que alcanza todo este material. 
La 9-4 muestra el dano mecAnico de la mezcla C-2 considerando que la porci6n de agregado 
p6treo susceptible de dano no colabora en la capacidad soporte de la mezcla, es decir se 
encuentra con un dano inicial del 100% 
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Sin dafio 60% dano 70% dano 

80% dano 90% dafio 100%dafto 



Figura 5: Curvas de respuesta del ensayo Marshall para la mezcla C-2 0, 6 < d h < 1 



Ensayo Marshall C-l 




Sin dafio io% dafto 20% dafio 

30% dafto 40% dafto 50% dafio 



Figura 6: Curvas de respuesta del ensayo Marshall para la mezcla C-l 0 < d h < 0, 5 
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Ensayo Marshall C-l 




Sin dafio 60% dafio 70% dafio 

80% dafto 90% dafio 100% dafio 



Figura 7: Curvas de respuesta del ensayo Marshall para la mezcla C-l 0, 6 < d h < 1 



C-l y C-2 - Carga Mdxima 




Figura 8: Evoluci6n de la capacidad portante en funcidn del dafio higrom6trico, para un 
desplazamiento prescrito de 1,8mm. 



l'n liiodelo constitutive) olastoplastico acoplado con dano mecdnico e higrometrico 
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7 CONCLUSIONES 

En el presente trabajo se muestra la teoria de mezclas sometida a efectos mec&nicos e 
higrometricos para el andlisis de materiales compuestos sometidos a fendmenos de dano y 
plasticidad. Ademas, tambien se presenta una correccidn a la teoria de mezclas cl&sicas, que 
permite tratar materiales de matriz compuesta. Por ultimo, se aplica esta teoria a mezclas 
bituminosas, utilizadas para la construction de pavimentos de carreteras. En este ultimo 
caso se ofrece una comparacidn con el ensayo Marshall obtenido en laboratorio, donde se 
aprecia una buena concordancia entre los resultados numericos con los obtenidos a partir 
de ensayos. 
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RESUMEN 

En este articulo se describe un modelo de analisis no lineal paso a paso en el tiempo, de 
estructuras de hormigon construidas evolutivamente, en condiciones de servicio y en situacion de 
estado lfmite ultimo. La consideration del proceso constructive evolutivo permite extender el 
modelo para simular intervenciones tales como demoliciones parciales, reparacion mediante 
inyeccion de fisuras o refuerzo mediante encolado de chapa o pretensado exterior. El analisis 
incremental permite evaluar la capacidad resistente residual y el comportamiento en servicio de la 
estructura reparada, siendo de gran ayuda para la selection del sistema mas eficaz en terminos 
economicos y estructurales. Se presenta, a modo de ejemplo, los analisis realizados para estudiar la 
necesidad de reparacion del puente de Las Pilas en Huesca, simulando el estado de fisuracion y 
deformaciones del mismo y los efectos de la introduction de un pretensado exterior. 



ABSTRACT 

An analytical model for the nonlinear and time-dependent analysis of segmentally erected 
concrete structures under service and ultimate load conditions is presented. The model allows to 
simulate the effects of repair and strengthening interventions such as cracks grouting, partial 
demolition, bonding of steel plates or introduction of external prestressing, among others. The 
incremental analysis procedure allows to evaluate the residual ultimate load and the service 
behavior of the retrofitted structure, being useful to select the most suitable retrofit system in terms 
of economy and structural efficiency. The analyses carried out to simulate the effects of introducing 
external prestressing on the cracking and deformation state of Las Pilas bridge in Huesca (Spain), 
are shortly presented as a matter of example. 
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1 INTRODUCTION 

La intervention en estructuras existentes de edification y de obra publica es un hecho cada vez 
mas frecuente. Por una parte, la agresividad ambiental, la falta de planteamientos sobre durabilidad 
en proyecto o el escaso mantenimiento de las estructuras hacen necesario reparar los danos 
producidos por el deterioro de los materiales a lo largo de la vida util de la estructura, a fin de evitar 
su progresion. Por otra parte, la aparicion de acciones imprevistas o cualquier modification en 
cuanto a geometria, cargas o forma de apoyo, orientada a adaptar la estructura a nuevas condiciones 
de uso, suele requerir la ejecucion de un refuerzo de la misma. La figura 1 muestra el estado de una 
pila hueca del puente que cruza el embalse de La Baells, como consecuencia de una presion 
hidrostatica excesiva, la cual hubo de ser reparada y reforzada. 




Figura 1. Estado de una pila del viaducto de La Baells. 



En la actualidad se han desarrollado numerosas tecnicas de reparation y refuerzo de estructuras 
de hormigon armado y pretensado, las cuales quedan sintetizadas en las referencias (FIP, 1994; 
Roca et al 1995; Casas et al. 1996; CEB, 1983; Allen et al 1987; Priestley et al 1996). Sin 
embargo, para que un refuerzo sea teen icamente eficaz y economicamente aceptable, debe partir de 
un buen conocimiento del estado de la estructura que se refuerza y, sobre todo, debe estudiar los 
efectos que la intervention introduce sobre su comportamiento en servicio y sobre su capacidad 
residual. Elio implica enormes dificultades debido al desconocimiento de multitud de variables, 
como la forma de construction, la historia de cargas y otros aspectos ligados a la vida util de la 
estructura. Por ello suele acudirse a calculos simplificados, generalmente basados en hipotesis de 
linealidad, que pueden conducir a sobredimensionar el refuerzo. Pero ademas, muchas veces se 
producen grandes vacios de conocimiento sobre la eficacia real de la intervention. 

La capacidad de los ordenadores actuales y el desarrollo de los metodos numericos permiten 
disponer de potentes herramientas de analisis estructural capaces de simular con cierto realismo el 
comportamiento instantaneo y diferido no lineal por los materiales y la geometna, y los efectos del 
proceso construed vo evolutivo en esquema longitudinal y section transversal. Estos metodos 
quedan sintetizados en las referencias (Meyer et al 1985; Bazant, 1982; Mari, 1984; Abbas, 1993; 
Cruz et al. 1998). 

Sin embargo se echa en falta la disponibilidad de modelos generales de analisis capaces de 
reproducir tanto el proceso de deterioro de las estructuras como las actuaciones de reparation o 
refuerzo. Entre los fenomenos asociados al deterioro estan la fisuracion del hormigon, la 
degradation de las propiedades tenso-deformacionales de los materiales debido a procesos ffsico- 
qufmicos (corrosion, aumento de la porosidad, expansividad, microfisuracion, acciones ciclicas), la 
perdida de adherencia entre materiales y la perdida de section, entre otros. 
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Las actuaciones de reparacion o refuerzo pueden ser tan variadas como la inyeccion de fisuras 
con resinas o morteros, la aplicacion de pretensado exterior, el recrecido con hormigon y armaduras, 
el encolado de chapa metalica, el zunchado de pilares mediante presillas o el refuerzo con materiales 
compuestos. Todo ello con la posibilidad de apear o descargar la estructura, modificar su esquema 
estatico, introducir lastres u otras operaciones de diversa indole. 

En este articulo se describe un modelo general de analisis no lineal paso a paso en el tiempo, 
actualmente en desarrollo en el Departamento de Ingenieria de la Construction de la UPC, que 
permite simular algunas de las operaciones anteriores. Se presenta, a modo de ejemplo de aplicacion 
del modelo, el estudio realizado sobre un puente real (el puente de las Pilas en Barbastro, Huesca). 
En tal caso, el modelo permitio comparar la respuesta estructural frente a diversas soluciones de 
reluerzo, y fue determinante en la decision sobre la solution adoptada. 

2 CARACTERISTICAS DEL MODELO DE ANALISIS. 



2.1 Generalidades y campo de aplicacion. 

Se trata de un modelo de analisis de porticos pianos o espaciales basado en una idealization de 
las estructuras mediante elementos lineales tipo barra con 6 grados de libertad por nodo, 
discretizando la seccion transversal en filamentos. 

Permite considerar de manera realista las propiedades tenso-deformacionales de los materiales y 
su influencia en la respuesta estructural, gracias al analisis no lineal. Asf, es posible reproducer la 
formation de fisuras y estudiar su influencia en el comportamiento en servicio y bajo niveles 
avanzados de carga. Tiene en cuenta tambien los efectos de segundo orden, de notable influencia en 
el comportamiento de estructuras esbeltas. 

El modelo permite considerar la evolucion de las propiedades del hormigon en el tiempo, los 
efectos estructurales de las deformaciones diferidas de retraction y fluencia del hormigon y de la 
relajacion del acero de pretensado y la interaction con el acero o con otros hormigones que tengan 
dilerentes edades de fabrication y puesta en carga. Las redistribuciones de esfuerzos a nivel 
estructura, de tensiones a nivel seccion y el incremento diferido de flechas son, entre otros 
resultados del analisis. 

Permite considerar la evolucion de la configuration estructural durante el proceso constructivo, 
por ejemplo las variaciones en esquema longitudinal, el hormigonado de las secciones en diferentes 
etapas, la simulation de rotulas que pueden ser bloqueadas durante el proceso de construccion, la 
colocation y elimination de cimbras, tirantes u otro tipo de apoyo provisional, el destesado y 
retesado de tendones y tirantes, etc. 

De manera general, mediante la utilizacion del modelo se pueden estudiar el comportamiento en 
servicio y la evolucion hasta rotura de porticos pianos y espaciales, bajo cargas instantaneas o 
duraderas, la distribucion de reacciones, esfuerzos, tensiones y deformaciones durante la 
construccion, los efectos del proceso constructivo en el comportamiento en servicio y en la 
capacidad portante y, como se explicara mas adelante, la respuesta estructural frente a determinados 
tipos de intervenciones. 

HI programa que permite la utilizacion practica del modelo matematico se denomina CONS, y ha 
sido desarrollado por A. Mari a partir de no previamente desarrollado para analisis no lineal, 
llamado PCH3D (Abbas & Scordelis, 1993). Esta confeccionado en forma modular, con lo cual 
permite su continua evolucion, introduciendo nuevos fenomenos mediante modificaciones parciales 
que lo ajustan cada vez mas al comportamiento real de las estructuras. 

2.2 Cargas exteriores 

Todas las cargas se aplican sobre los nodos de la estructura en instantes precisos, conocidos 
como escalones de carga, y esta puede dividirse al efecto en varios incrementos. Para la realizacion 
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del analisis paso a paso en el tiempo se considera que las cargas se aplican solo en los escalones de 
tiempo, por lo cual la tension en cada punto de la estructura es constante durante cada intervalo. 

El pretensado se introduce como un sistema autoequilibrado de cargas y momentos aplicados en 
los nodos, obtenido por equilibrio del cable. Aunque el trazado de los tendones puede ser rectilfneo, 
quebrado o parabolico, este se discretiza en segmentos rectilineos que cubren un solo elemento y 
concentran la interaccion con el hormigon en los puntos de quiebro. El modelo permite la aplicacion 
de la fuerza de pretensado desde un extremo del tendon o desde los dos. El programa evalua la 
fuerza a lo largo del tendon teniendo en cuenta las perdidas instantaneas de tension por rozamiento y 
penetracion de curias y, por equilibrio del cable genera un sistema de cargas nodales. Las perdidas 
por acortamiento elastico se obtienen como resultado del calculo, en funcion de la deformacion de 
la estructura. Si se trata de armaduras postesas, el tendon que se tesa no se deforma con la estructura 
hasta que se ha anclado. En el caso de armaduras pretesas la adherencia establece la compatibilidad 
de deformaciones, dando lugar a perdidas por acortamiento elastico al transferir. 

El programa admite desplazamientos impuestos en los apoyos, asf como deformaciones 
impuestas debidas a la accion termica ambiental o a la retraccion del hormigon. 

2.3 Esquema estructural. 

La estructura se discretiza longitudinalmente mediante elementos finitos Hermitianos tipo barra 
de 3 nodos, con seccion constante y forma arbitraria. El elemento tiene 13 grados de libertad, 12 en 
los extremos y uno adicional de deformacion axial en el centro (que se elimina por condensacion 
estatica a nivel del propio elemento), necesario para obtener correctamente la matriz de rigidez a 
flexion y tener en cuenta la variation del eje neutro por fisuracion. 

Las condiciones de contomo se reproducen mediante apoyos, definidos con 6 resortes cuya 
rigidez define el usuario y puede variar a lo largo del tiempo. Dichos resortes controlan los giros y 
los desplazamientos de la estructura en estos puntos. Estas condiciones pueden variar en el tiempo, 
permitiendo reproducir los procesos de cimbrado, descimbrado, colocacion de nuevos apoyos y 
eliminacion de los existentes. 

La colocacion de nuevos apoyos se consigue restringiendo el desplazamiento de la estructura en 
ese sentido para todas las cargas que actuen a partir de ese momento. La eliminacion de una 
condition de apoyo existente consiste en aplicar la carga en el apoyo eliminado, como una fuerza 
desequilibrada igual y contraria que actua en el siguiente escalon de carga o de tiempo sobre la 
nueva configuration estructural. 

La liberation de grados de libertad extremos de los elementos, permite considerar la presencia de 
rotulas entre elementos, el deslizamiento relativo entre ellos, etc. Estos vfnculos intemos pueden 
variar a lo largo del proceso constructive o durante la evolution de la estructura. 

2.4 Comportamiento a nivel seccion 

Se adopta la hipotesis deformacion plana de la seccion, que puede expresarse asf: 



A £(y, z ) = Ae m - z • A c v + y ■ Ac z 



( 1 ) 



Donde: 

Ae (y,z) incremento de deformacion en el filamento situado en las coordenadas (y, z) 
respecto al sistema de referencia de la seccion 
Ae m incremento de deformacion en el centro de referencia de la seccion 

y,z coordenadas del filamento (medidas desde el centro de referencia) 

Ac y , Ac z incrementos de curvatura 

Se considera adherencia perfecta entre los materiales que forman la seccion. Esto se cumple 
tambien para el acero de pretensado una vez se considera adherente. 

Se supone que cada elemento es un prisma cuya seccion transversal de hormigon es un solido 
discretizado en filamentos. Las caracterfsticas geometricas de la seccion neta resulta de la 
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integracion de las areas de los filamentos que la componen y cuando alguno de estos filamentos no 
existe (ya sea por no haberse construido aun, haber fisurado, etc.) la contribucion de su area 
diferencial es cero. Esto permite tener en cuenta la aparicion o desaparicion de filamentos segun lo 
exija el proceso constructive y considerar para cada uno de ellos la ecuacion constitutiva que le 
corresponde. 

2.5 Comportamiento a nivel punto. 

La delormacion uniaxial total de un punto en un instante t n puede descomponerse en: 



n total . nm 

£ n ~ £ n + e n 



( 2 ) 



donde: 

£ n Es la deformacion mecanica, entendida como tal, la deformacion instantanea debida a una 
tension aplicada. 

£ n Es la deformacion no mecanica, entendida como tal toda deformacion no correspondiente a 
una tension instantanea. En ella se encuentran englobadas las deformaciones debidas a 
fluencia, retraccion, variaciones termicas y, de forma ficticia, el envejecimiento del 
hormigon. En el caso del acero esta se reduce a la deformacion termica. 

En base a ello, el calculo de las tensiones y deformaciones en un escalon de tiempo t n se obtiene 
mediante el procedimiento siguiente: 

1. Se calcula la deformacion total en el instante t n acumulando la deformacion obtenida en el 
iotervalo t n -t n .j , sobre la total del escalon anterior: 



e n =e n _ l +Ae n ( 3 ) 

2. Se obtiene el incremento de deformacion no mecanica entre los escalones t n .\ y t n como: 

c =< + <+< + £ (4) 



3. La deformacion no mecanica total en el instante t n se obtiene acumulando el incremento en este 
intervalo sobre el total anterior. 



£ 



nm 

n 



-£ 



nm 
n - 1 



+ Ae 



nm 

n 



4. La deformacion mecanica se obtiene entonces deduciendo de la total la no mecanica. 



,m total nm 

n ~ k/t 



(5) 



( 6 ) 



Las tensiones se determman a partir de la deformacion mecanica como: 

^=/(C) = /(C d/ -C) 



(7) 



Donde: 

£ n Deformacion uniaxial mecanica de un punto en el instante t n 

£ n Deformacion uniaxial no mecanica de un punto en el instante t n 

£ n -i Deformacion total en el instante t n .j 

A£ n Incremento de deformacion total en el intervalo entre t n _i y t n 

£c Incremento de deformacion por fluencia del hormigon entre tn_j y t n 

£c r Incremento de deformacion por retraccion del hormigon entre t^, y t n 
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ej Incremento de deformacion por variacion de temperature entre t n .i y t n 

£c a Incremento de deformacion por envejecimiento del hormigon entre t n .i y t n 

£st Incremento de deformacion en el acero por variacion de temperature entre Vi y t n 

£n-i nrn Deformacion no mecanica acumulada hasta el instante t„_i 

A£nnm Incremento de deformacion no mecanica en el intervalo entre t„.i y t n 

Para llevar a cabo este analisis paso a paso, se considera el nivel de tensiones constantes durante 
un intervalo de tiempo y posteriormente se obtienen los nuevos incrementos por superposicion al 
variar la historia de cargas. Para historias de carga decrecientes, el principio de superposicion 
subestima la recuperacion de fluencia. 

2.6 Comportamiento instantaneo de los materiales 

El hormigon se considera sometido a un estado uniaxial de tensiones normales. La ecuacion 
constitutiva adoptada para el hormigon es parabolica y lineal y comprende 1 1 estados diferentes, 
como indica la figure 2. En ellos se determina la tension del hormigon a partir de la deformacion 
mecanica y el modulo de deformacion del material en el instante que se analiza. Se considera que ha 
fisurado si la tension de traccion excede su resistencia a traccion en ese instante y plastifica cuando 
su deformacion mecanica excede la correspondiente a su resistencia maxima a compresion, 
comenzando una rama de reblandecimiento. 



ct 




Figura 2. Ecuacion constitutiva adoptada para el hormigon. 

La fisuracidn se considera distribuida ("smeared cracking"), por Io que la contribucion del 
hormigon traccionado entre fisuras se tiene en cuenta a traves de una rama de descarga hiperbolica 
en el diagrama o-£ del hormigon en traccion. 

El acero de las armaduras pasivas tiene un comportamiento reologico muy diferente al hormigon, 
no existiendo retraccion y no siendo apreciables los efectos de fluencia o relajacion. Por ello se 
supone su comportamiento constante en el tiempo, basandose en un diagrama bi lineal simetrico 
respecto al origen con descarga, como el que muestra la figura 3. 

El comportamiento del acero de pretensado se Simula a partir de un grafico multi lineal, 
suponiendo siempre que la descarga y recarga ocurren con el modulo de deformacion inicial. En la 
figure 4 se representa el comportamiento de los tendones de acero pretensado. 
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Figura 3. Ecuacion constitutiva adoptada para el acero pasivo. 



Tension [MPa] 




Figura 4. Ecuacion constitutiva del acero de pretensado. 



2.7 Comportamiento diferido de los materiales 

La deformacion de fluencia (deformaciones tensionales diferidas del hormigon) se tiene en 
cuenta considerando las hipotesis de viscoelasticidad lineal, o sea, linealidad entre tensiones y 
deformaciones diferidas y el principio de superposicion de Boltzman que se puede expresar como: 

£« (0 = <H t 0 )c(t,t 0 ) + Y i=i Act (t , ) • c(t,t, ) (g) 



Donde: 

£cc(t) Deformacion de fluencia 

a(to) Tension aplicada en el instante to 

Ao(tj) Increments de tension aplicados entre los instante to y t 

c(t,ti) Valor de la fluencia especi'fica en el instante t para una carga aplicada en t, 
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Para evitar la necesidad de almacenar toda la historia de cargas, se utiliza para esta funcion de 
fluencia un desarrollo en forma de serie de Dirichlet: 

c«.<- r) = i;„«,(r)|l- e - WK -'’| (» 



Donde: 

aj (x) son funciones dependientes del instante de carga, T, cuyo valor influye en el 

coeficiente final de fluencia 

Ai son constantes que gobieman el desarrollo temporal de la fluencia 

Estos parametros se obtienen automaticamente en el programa mediante un ajuste por mfnimos 
cuadrados a la curva experimental o analftica que se desea reproducir. 

De esta manera, utilizando el metodo paso a paso en combinacion con la serie de Dirichlet, se 
logra que cada nuevo incremento de deformacion de fluencia dependa solamente de la tension y de 
unas variables de estado correspondientes al ultimo incremento de tiempo y no de toda la historia 
tensional. 

Las variaciones de tension del acero de pretensado producidas por la relajacion bajo una 
deformacion constante, se tienen en cuenta conjuntamente con el resto de deformaciones producidas 
en los escalones de carga o de tiempo por compatibilidad con el hormigon de esa fibra mediante el 
metodo de la tension inicial ficticia, obtenida mediante la siguiente ecuacion basada en datos 
experimentales: 



A = i_ logf 

fu 



10 



Ll 

fy 



-0,55 



> 0.55 

fy 



(10) 



Donde: 

f si Tension inicial ficticia. Aquella con la que debfa tesarse inicialmente el tendon para que su 
tension en el tiempo actual t fuese f s 

f s Tension actual en el acero de pretensado. Es la diferencia entre la tension inicial y las 
perdidas (por relajacion y por deformaciones elasticas en los escalones previos) 

f y Tension maxima del acero de pretensado 

2.8 Estrategia de analisis no-lineal en el tiempo 

El proceso de calculo se realiza mediante un proceso paso a paso que implica los cuatro bucles 
que se explican a continuacion, ordenados de forma que cada uno de ellos encierra a los siguientes: 

Etapas constructivas. Se definen como tales todas aquellas situaciones en las que se produce 
cualquier modificacion de la geometrfa, las cargas o las condiciones de vinculacion. 

Pasos de tiempo. El tiempo transcurrido de una etapa constructiva a otra se divide en intervalos, 
separados por escalones, en cada uno de los cuales se realiza un analisis frente a las acciones 
procedentes de las deformaciones diferidas, actual izando las propiedades de los materiales, si varfan 
en el tiempo. 

Aplicando el principio de los trabajos virtuales, las deformaciones no mecanicas se transforman 
en un sistema de fuerzas exteriores ficticias equivalentes que producirfan el mismo estado de 
deformaciones en el elemento. El vector de fuerzas no mecanicas producido en el intervalo de 
tiempo t n .j a t n por las deformaciones no mecanicas es: 

AR” m =\b t E t - Ae”"‘dv ( 11 ) 

El sistema de fuerzas exteriores actuantes sobre la estructura en el instante t n sera entonces: 
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AR n =AR' n+ AR:" + ARl { 



( 12 ) 



donde : 

AR n Vector de fuerzas actuantes al final del intervalo que va del instante t^ al t n 

Ae n mn Deformacion no mecanica producida en el intervalo t n .i al t n 

AR n e Vector de cargas exteriores que aparece en el intervalo t n .i al t n 

ARn Vector de cargas debido a las deformaciones no mecanicas ocurridas en el intervalo 

ARnu Vector de fuerzas residuales o desequilibradas, que se explica mas adelante. 



Los filamentos que componen la section transversal se mantienen a una tension constante 
durante cada intervalo de tiempo. El analisis parte de que se conoce el estado de tensiones, 
deformaciones y desplazamientos del escalon previo, y se determinan los incrementos en el actual. 

Escaiones de carga. En cada paso de tiempo, se realiza un proceso incremental en el cual la 
carga exterior total se divide en escaiones. Elio permite obtener la respuesta estructural bajo cargas 
crecientes. Este proceso es indispensable para conocer de manera realista la respuesta estructural, 
pues esta depende de la historia de cargas, debido al comportamiento no lineal del hormigon y del 
acero. Tambien existe la option de realizar un proceso incremental con desplazamiento controlado, 
lo que puede permitir reproducir fenomenos de inestabilidad o reblandecimiento. 

Iteraciones de equilibrio. En cada uno de los escaiones de carga se realizan iteraciones hasta 
conseguir la satisfaction de las condiciones de equilibrio. El proceso iterativo se realiza por los 
metodos de Newton-Raphson o Newton-Raphson modificado. En cada iteracion se llevan a cabo las 
siguientes operaciones: 



1 . 



A partir de la geometna y las propiedades de los materiales actualizadas se forma la matriz de 
rigidez tangente para cada elemento y se ensambla la de toda la estructura con la matriz de 
transformation actualizada. 

Se resuelve el sistema incremental AR= K r Ar, se obtienen los incrementos de desplazamiento 
Ar y con ellos los incrementos de deformacion de cada filamento Ae=B Ar. Estos Ae se ahaden a 
los previos totales para obtener las deformaciones totales actualizadas del hormigon, acero 
pasivo y de pretensado. 

Se acumulan los incrementos de desplazamientos Ar, obteniendo los desplazamientos totales de 
los nodos. Se actualiza la geometna de los elementos y su matriz de deformacion. 

A partir de las ecuaciones constitutivas de los materiales se obtienen las tensiones. 

Se integran estas tensiones obteniendose a traves del principio de los trabajos virtuales, el vector 
de cargas sobre los nodos, energeticamente equivalente al campo de tensiones. El ensamblaje de 
este vector para toda la estructura permite obtener el vector respuesta R 1 . 

Se obtiene el vector de fuerzas desequilibradas R u , como la diferencia entre el vector de cargas 
extemas totales R c y el vector respuesta, es decir R u = R e -R' . 

Se iguala AR = R u y se comienza una nueva iteracion desde 1, hasta que ocurra la convergence 
del proceso, bien porque las fuerzas desequilibradas (residuos) no superen un cierto valor 
establecido, bien porque las variaciones de energfa total del sistema sean suficientemente 
reducidas. 

Una vez Iograda la convergencia se introduce un nuevo incremento de carga y se repite todo el 
proceso iterativo. Cuando se han introducido todos los escaiones de carga dei instante considerado 
se procede al siguiente paso de tiempo en que se divide el intervalo de tiempo entre dos etapas 
consecutivas. El proceso siempre es acumulativo en resultados, de forma c „e se puede reproducir la 
historia real de construction y de cargas de la estructura. 

Dado que se trata de una herramienta de analisis no lineal, requiere partir del conocimiento de 
datos geometricos, materiales, cargas y tambien de las armaduras activas y pasivas. 



2 . 
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3 TRATAMIENTO DE LOS PROCESOS DE REPARACION Y REFUERZO 

La reparacion y el refuerzo de una estructura son intervenciones en la misma que suelen xmplicar 
modificaciones en su geometri'a, en las cargas, en el comportamiento de los materiales o en las 
vinculaciones. El planteamiento paso a paso del modelo matematico desarrollado permite tener en 
cuenta estas intervenciones. Para ello se especifican algunos aspectos en la lectura de datos, como 
s° n: 

- Cada filamento de una seccion transversal puede ser de un hormigon diferente, del cual se 
define el dfa en que se fabrica, el dfa en que se demuele y el dfa en que se repara inyectando las 
fisuras. 

- Para cada aimadura pasiva se define el dfa en que se coloca en las estructuras y el dfa que se 
demuele o retira. 

- Cada tendon de pretensado puede tesarse, destesarse o retesarse a lo largo del tiempo. Por tanto, 
en cada etapa constructiva se definen aquellos tendones que se tesan, retesan o destesan. 

- De cada elemento se especifica que dfa se ensambla en la estructura y que dfa se retira. 

La colocation de una nueva armadura, filamento de hormigon o barra estructural implica 
unicamente incorporar su rigidez a la estructura y activar algunos grados de libertad. 

La retirada o demolition de un elemento estructural requiere la introduction, como cargas, de los 
esfuerzos en extremos de la barra que se retira, y la desactivacion de algunos grados de libertad. 

A partir de estos principios, se pueden simular, entre otros los siguientes procesos de 
intervention: 

Saneado de recubrimientos. Suele implicar no solo una perdida de seccion de hormigon sino 
tambien de acero de las armaduras pasivas. Su simulation consiste en especificar que el hormigon 
del recubrimiento se demuele a una edad determinada. La perdida de seccion se Simula eliminando 
algunas armaduras (aunque estrictamente no sea asf). 

Inyeccion de fisuras. Si un hormigon fisurado se repara con esta tecnica se supone que vuelve a 
resistir tracciones. Por tanto se trata de una modification de la ecuacion constitutiva. 

Recrecido con mortero y armaduras o mediante encolado de chapa. Puede simularse 
especificando que un conjunto de filamentos de hormigon y armaduras se colocan a una 
determinada edad. La chapa se considera una armadura mas, con mayor superficie de contacto, de 
un acero de calidad distinta a las barras corrugadas. 

Pretensado exterior Se Simula especificando que ciertos tendones se introducen y tensan en un 
instante determinado. El pretensado introduce un sistema de cargas autoequilibrado que puede 
cerrar fisuras y modificar el estado de tensiones, deformaciones y reacciones de la estructura. 

Apeos, desapeos y movimientos impuestos, Se simulan gracias a la posibilidad de modificar a 
lo largo del proceso paso a paso, las condiciones de vinculacion de la estructura. 



4 EJEMPLO: EL PUENTE DE LAS PILAS EN BARBASTRO (ESPANA) 



4.1 Introduction. 

El puente de “Las Pilas” sobre el no Cinca se encuentra ubicado cerca de Barbastro, en la 
provincia de Huesca, (Espafia). Su esquema es el de una viga continua de seccion cajon de 5 vanos, 
siendo los dos vanos extremos de 30 m. de luz y los 3 centrales de 60 m, como indica la fotografia 
de la figure 5. El canto es variable, de 1.60 m en centra de vano y estribos, y 2,40 m sobre apoyos 
intermedios. 
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Figura 5. Vista general del puente de Las Pilas, en Barbastro (Huesca). 

La particularidad del puente es que se trata de un tablero continuo en seccion cajon constituido 
por vigas prefabricadas y losa superior hormigonada in situ . La seccion de las vigas es en artesa, 
existiendo diafragmas en los extremos de las mismas y sobre apoyos. En total existen tres tipos de 
piezas: los ’’sombreros", o piezas de 20 m de longitud que se situan sobre las dos pilas centrales, las 
vigas extremas de 30 m de vano y 10 m de voladizo, y las tres vigas centrales, de 40 m de longitud, 
que se apoyan sobre los sombreros y/o sobre las vigas extremas (Ver figuras 6 y 7). La continuidad 
se consigue gracias al postesado de las juntas, previamente inyectadas con mortero de resina, y a la 
losa superior de hormigon armado, ejecutada in situ mediante carro de avance. 




El proceso constructive resulta enormemente evolutivo, tanto en esquema longitudinal como 
transversal. Por una parte se disponen apeos provisionales en los sombreros durante el montaje. Por 
otra parte se hormigona la losa por etapas: primero la zona de apoyos, para dar suficiente resistencia 
y rigidez a la seccion y despues la zona de centro de vano. Ademas transversalmente se hormigona 
en dos etapas: primero entre almas, para cerrar el cajon, y luego los voladizos. Por ultimo, se 
proporciona continuidad estructural, modificando, por tanto las vinculaciones entre elementos. 

Este proceso tan evolutivo, junto a las deformaciones de fluencia y retraction del hormigon, dan 
lugar a notables flechas diferidas y a redistribuciones de esfuerzos longitudinales y de tensiones 
entre viga y losa, hormigon y acero, a lo largo del tiempo, que deben ser consideradas para conocer 
la respuesta estructural del puente. 
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AJ poco tiempo de finalizar la construccion, se midieron flechas del orden de 14 cm en los tres 
vanos centrales, valor que superaba ligeramente los valores lfmites aceptables (1/500 bajo cargas 
permanentes, es decir 12 cm), y resultaban apreciables a simple vista e incluso generaban una cierta 
incomodidad al conductor al circular sobre el puente. Elio provoco la realization de un estudio 
detallado del puente para conocer su estado tenso-deformacional y los efectos de posibles 
intervenciones sobre el mismo. 
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Figura 7. Esquema longitudinal. 



4.2 Alcance del analisis realizado. 

Una inspection in situ permitio observar el estado de fisuracion de las vigas. En las secciones de 
apoyo se apreciaban fisuras por flexion en la zona de moment os negativos de hasta 1,45 m de 
profundidad maxima a 2.5 m del eje del apoyo. Esta fisuracion se comenzaba a inclinar al alejarse 
del apoyo por efecto del esfiierzo cortante. La fisuracion desaparecia aparentemente una vez 
finalizado el voladizo, o sea en las vigas intermedias a 10 m de los apoyos. 

Antes de valorar las posibles actuaciones a llevar a cabo, se imponfa identificar las causas de las 
flechas, verificar su estado actual e intentar reproducirlo lo mas precisamente posible. Para ello se 
simulo, mediante el modelo descrito anteriormente, el proceso construed vo y de puesta en servicio 
del puente, contrastando los resultados obtenidos, al nivel de flechas y fisuracion con los que se 
habian observado y medido en la estructura construida, obteniendo resultados bastante 
concordantes. 

Las conclusiones extraidas en cuanto a fisuracion y flechas fueron las siguientes: La fisuracion 
aparecida sobre los "sombreros", en la zona de momentos negativos, estaba prevista en el proyecto, 
y era debida a la ausencia de pretensado en dicha zona. 

Sin embargo, las flechas observadas eran superiores a las previstas en proyecto. Las causas de su 
aparicion se imputaron a: la escasa colaboracion de la losa superior por efectos de la prefisuracion 
del hormigon por retraction, fruto de la coaccion ejercida por la viga a la libre deformation de la 
losa, a la curvatura impuesta por el pretensado de continuidad excentrico y a la enorme diminution 
de rigidez de la zona de apoyos que supuso la fisuracion. 

Una vez reproducido el estado real del puente se plantearon dos tipos de intervention: 
a) La disposicion de un pretensado exterior, dispuesto en el interior del cajon, que contribuyera a 
la recuperation de las flechas durante un perfodo de 20 anos y el sellado de fisuras. Se 
estudiaron las tres variantes que se esquematizan en la figura 8, consistentes, basicamente en: 
al) Disposicion de cuatro tendones de 24 <() 0,6", con trazado recto sobre la zona de momentos 
negativos en los sombreros, con una fuerza total de pretensado de 18000 kN. 
a2) Disposicion de tres tendones de 24 § 0,6”, con trazado poligonal con quiebros a cuartos de 
la luz, con una fuerza total de pretensado de 13500 kN. 
a3) Disposicion de tres tendones de 24 <|> 0,6", con trazado poligonal en los centros de vano, 
con una fuerza total de pretensado de 13500 kN. 
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ANCLAJE 4 torones de 24 0 0.6" 




SOLUCI 6 NA .1 




SOLUCI 6 NA .2 




SOLUCI6NA.3 



Figura 8. Soluciones de refuerzo con pretensado exterior. 
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b) Una escarificacion del pavimento, disposition de hormigon ligero, repavimentacion de la 
calzada, sellado de juntas y obras de menor envergadura en barandillas e imposta, que 
mejoraran sensiblemente el estado de la rasante y el aspecto visual del puente. 

Evidentemente la diferencia de coste entre las soluciones con pretensado exterior y la solucion b) 
era enorme, ya que las primeras implicaban realizar aberturas en la losa superior, disponer 
desviadores, hormigonar macizos de anclaje, disponer y tesar tendones en el interior del cajon y 
operaciones de reparation y acabado. 

A fin de disponer de information suficiente para tomar una decision se decidio, aparte de tomar 
otras medidas y realizar calculos mas simplificados, simular, mediante el modelo descrito, la 
respuesta estructural del puente en condiciones de servicio, a corto y largo plazo, en cada uno de los 
casos de intervencion estudiados. Asimismo se decidio simular un proceso de sobrecarga hasta 
rotura para conocer la capacidad resistente a flexion. 

4.3 Simulation numerica. 

A partir de los datos disponibles sobre la fabricacion de las vigas, fuerzas e instantes de tesado y 
de la description del proceso constructivo real, se ajustaron las fechas de las etapas constructivas y 
pasos de tiempo a suministrar al modelo, procediendose a simular el proceso constructivo de forma 
precisa. Se analizo solo medio puente por simetria, discretizando la estructura en 60 elementos de 2 
m, cada uno de ellos con section de hormigon y/o armadura constante. 

La caracterizacion de sus propiedades mecanicas se hizo segun las especificaciones del CEB. Se 
utilizaron hormigones de dos tipos diferentes: para la losa un hormigon de 30 MPa de resistencia 
caracteristica a los 28 dfas, y para la viga uno de 50 MPa. Cada hormigon se considera con edad 
cero el dia de su fabricacion. La armadura pasiva era de acero corrugado de 500 MPa de 1 unite 
elastico. La armadura de pretensado es de acero superestabilizado 138-79 segun Euronorma, con un 
lfmite elastico de 1700 MPa y una relajacion del 2.5% a las 1000 horas. 

Aparte del peso de la viga y de la losa se consideraron las cargas debidas a los carros de 
hormigonado, cargas muertas estandar de trafico ( 4 kN/m 2 + carro de 600 kN) y acciones debidas a 
repavimentado (2,4 kN/m 2 por cada centimetro de espesor) y colocation de hormigon ligero (1,5 
T/m 3 ), dispuesto con grueso variable parabolicamente de valor maximo en los centros de vanos. Se 
considero un espesor maximo de hormigon ligero de 10 cm. Los tendones de pretensado se tensaban 
al 80 % de su tension de rotura. 

4.4 Resultados del calculo. 

La tabla 1 adjunta describe las flechas obtenidas en los calculos para las tres soluciones a base de 
pretensado exterior, para la solucion b), incluida la repavimentacion, y las obtenidas en caso de no 
actuar (inicial). Se observa que, entre las intervenciones con pretensado exterior, la solucion a-2 es 
la mas efectiva a efectos de flechas. Por ello, y por limitaciones de espacio, en lo que sigue se 
describiran con mas detalle unicamente los resultados de esta hipotetica actuation y de la solucion 
b). Tambien se indica en esta tabla la flecha de la rasante a largo plazo en la solucion, entendiendo 
como tal la resultante de restar a la flecha estructural el espesor del hormigon ligero dispuesto, cuyo 
peso se incluye en los calculos. Dado que en las soluciones a) no se dispone hormigon ligero, la 
flecha y la rasante coinciden. 



SOLUCION 


Flecha instantanea 


Flecha total largo plazo 


Rasante a largo plazo 


Vano 2 


Vano 3 


Vano 2 


Vano 3 


Vano 2 


Vano 3 


inicial 


133 


150 


168 


158 


168 


155 


a-1 


108 


132 


150 


142 


150 


142 


a-2 


95 


122 


130 


134 


130 


134 


a-3 


100 


128 


138 


141 


138 


141 


b 


147 


155 


180 


165 


80 


65 



Tabla 1. Flechas y estado final de la rasante en mm para diversas soluciones de intervencion. 
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Solucion a-2: pretensado exterior con trazado poligonal y quiebros a cuartos de la luz. 

Los efectos estructurales mas significativos, que pueden observarse en las figuras adjuntas son: 

1 . Flechas. Se recuperan de forma instantanea 3.8 cm en el vano 2 y 2.8 cm en el vano 3. A largo 
plazo se incrementa la flecha, respecto de la instantanea, en 3,5 cm en el vano 2, y 1,2 cm en el 
vano 3. La figure 9 muestra las deformadas obtenidas con esta solucion. 




Figura 9. Deformadas solucion a-2. 

2. Fisuracion. Tras pretensar se cierran fisuras en las vigas armadas en los voladizos y disminuye 
la profundidad de las fisuras. Tambien se cierran fisuras en la losa, situation que persiste a largo 
plazo. 

3. Momentos y tensiones- Aumentan ligeramente los momentos negativos (descontando los 
isostaticos) y disminuyen los positivos (aproximadamente en un 10%) al pretensar. Es decir el 
efecto de esta intervention es generar una ligera diminution de la reaction de los estribos y por 
tanto un aumento del momento negativo. No obstante debido a los efectos del pretensado en su 
zona directa de aplicacion, mejora el estado tensional. Las tensiones en el hormigon de la losa 
inferior aumentan ligeramente en las vigas 1 y 3, y tambien en las vigas 2 y 4, donde resulta 
favorable de care a fisuracion. Disminuyen las tensiones en las armaduras pasivas superiores de 
las vigas armadas (de 3300 a 2700 Kp/cm2 en la viga 3) 

Solucion b) repavimentacion. 

La solucion b) estudiada trata de minimizar la actuation sobre el puente, escarificando el 
pavimento y restituyendo la rasante prevista a largo plazo a base de hormigon ligero y 
repavimentando. Los resultados en terminos de flechas se muestran en la figura 10. La fisuracion 
resulta practicamente inalterada y, aunque la flecha aumenta ligeramente respecto de la solucion a), 
la rasante mejora sensiblemente. 
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Figura 10. Deformadas con la solucion b). 




Capacidad resistente a flexion 

Para cada uno de los cuatro casos estudiados se simulo un proceso de aumento de la sobrecarga 
de 4 kN/m 2 , hasta alcanzar la rotura por flexion. La capacidad de carga a flexion del puente, a pesar 
de la fisuracion existente, es muy alta en todos los casos, como lo demuestra el resultado del analisis 
no lineal que se muestra en la figura 11. En ella debe tenerse en cuenta que solo se aumenta la 
sobrecarga y el carro, pero que las cargas permanentes no estan mayoradas. No obstante, el factor de 
seguridad determinista a flexion del puente es superior al valor establecido por la normativa vigente, 
supuesto que la resistencia a cortante sea suficiente. 




- 0.10 - 0.15 -0 20 - 0.25 - 0.30 -035 -040 - 0.45 - 0.50 - 0.55 

Ftecte en centra vano 2 (metros) 



Figura 11. Curvas cargas-desplazamiento hasta rotura. 

Se observa, ademas, que respecto de la solucion b), el pretensado exterior aumenta en un 10% , 
aproximadamente, la capacidad resistente a flexion, excepto en la solucion a-1, donde incluso la 
reduce. Posiblemente ello sea debido a que la introduction del pretensado refuerza la zona de 
momentos negativos, pero solo afecta a los centros de vano, donde en ese caso se produce la rotura, 
a traves de efectos hiperestaticos del mismo signo que las cargas exteriores, y, por tanto, 
desfavorables. 

Como consecuencia de todo lo anterior se decidio adoptar la solucion b), es decir no realizar 
refuerzos con pretensado exterior, sino un sellado de fisuras, una escarificacion, disposition de 
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hormigon ligero y repavimentacion y otras actuaciones tendentes a mejorar el estado de la rasante, a 
reducir la percepcion visual de la deformacion, y mejorar la-comodidad del conductor. Igualmente 
se decidio realizar un seguimiento del estado de deformaciones y fisuracion del puente, tanto 
durante las obras de reparacion como a lo largo del tiempo. Las obras se realizaron con exito 
durante la primavera del 997. 



5 CONCLUSIONES Y PERSPECTIVAS 

• Se ha presentado un modelo general de analisis no lineal de estructuras de hormigon que 
permite considerar el comportamiento no lineal instantaneo y diferido de los materiales, el 
proceso constructivo y posibles actuaciones posteriores de refuerzo o reparacion. 

• A traves de un ejemplo real se ha mostrado la capacidad del modelo para simular diversos tipos 
de intervenciones y evaluar la eficacia estructural de cada una de ellas. Elio tanto en lo referente 
al comportamiento en servicio como a la capacidad resistente a flexion. Por este motivo resulta 
de gran ayuda para tomar decisiones respecto de la solucion idonea, permitiendo evitar 
intervenciones innecesarias y reduciendo su coste, en ocasiones de manera drastica. 

• FJ modelo se encuentra en proceso de constante mejora, siendo imprescindible introducir la 
interaccion-flexion cortante a fin de reproducir de forma realista determinados tipos de roturas 
en estructuras reforzadas, entre las cuales se encuentran las roturas por cortante-flexion, por 
perdida de anclaje, o por "pilling" de las platabandas encoladas. 

• Actual mente se utilizan materiales compuestos por resinas y fibras de carbono, como materiales 
de reparacion, dadas sus altas prestaciones. Se pretende a corto plazo incorporar estos 
materiales en el modelo, incluyendo no solo su comportamiento a nivel de ecuacion 
constitutiva, sino tambien su interaccion con el hormigon, es decir la adherencia imperfecta, a 
traves de la model izacion de la interfaz. 
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RESUMEN 



En este trabajo se desarrolla la estimacion de un modelo estocastico de la accion sfsmica, tomando 
como base la information disponible para la ciudad de Manizales (Colombia). El modelo utilizado 
es el de espectro evolutive (Yeh and Wen 1990). La caracterizacion de los parametros se realiza por 
me 10 de tecmcas de identificacion no lineal de parametros y metodos propios del proceso de sefiales. 

El modelo es usado luego para estimar la vulnerabilidad de los edificios sin diseno sismo-resistente 
de la ciudad. 



SUMMARY 



A evolutive stochastic seismic model, proposed by Yeh and Wen (1990), is calculated for the city 
anizales (Colombia) using the available seismic records. The calculation of the model parameters 
is performed by nonlinear identification and signal processing techniques. The model is then used for 

X buildings 6 Se ‘ SmiC fragiHty ° f a bUUding WUh n ° aS6iSmiC deSigU ’ aS re P resen tative of the city’s 



1. INTRODUCTION 

la f lngenlerl ' a sl 'smica ha habido desde sus comienzos un interes persistente por consid- 

debido a lal 8 aC -n°T C ° m ° la y espuestas estructurales desde una perspectiva probabilista. 
debido a las multiples causas de estocasticidad e incertidumbre presentes en las variables 

Z S °^ZlL d ZT Si b T‘ l0S “ d ' g “ dC CO “ tniCd “"' I” sobierna,, la pr toi ca 
tamhie y f g conceptos y parametros probabilistas de manera expb'cita o implfcita 
tambien es cierto que estos resultan abiertamente insuficientes para determinar de manera 
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clara la respuesta estructural en terminos probabilistas y, especialmente, la probabilidad de 
fallo. Para este proposito la ingenieria estructural dispone de varios metodos, los cuales 
pueden clasificarse en dos grandes grupos, a saber (cf. Casciati and Faravelli 1991): 

1. Metodos analiticos. 

2. Metodos seudo-estadisticos, tambien llamados metodos de Monte Carlo. 

Mientras los primeros enfocan el problema desde una perspectiva basada en una combi- 
nation de la teoria de probabilidades y la dinamica de estructuras, con lo cual se obtienen 
ecuaciones diferenciales de algunas medidas probabilistas, los segundos buscan generar una 
muestra artificial de las respuestas estructurales por medio de la solution de multiples prob- 
lernas deterministas con datos aleatorios correspondientes a la action sismica y a las variables 
estructurales estocasticas. En ambos casos se requiere de modelos probabilistas. Variables 
tales como la resistencia de los materiales, las dimensiones de los elementos, etc. pueden 
ser modeladas simplemente por medio de su funcion de distribution de probabilidad. En 
el caso de del movimiento del terreno se requiere un modelo estocastico de la variacion 
temporal de las ondas condicionada a la ocurrencia del fenomeno, ademas de la description 
probabilista de esta, que usualmente se asume de tipo Poisson. A partir de este modelo es 
posible generar los multiples acelerogramas artificiales requeridos por el metodo de Monte 
Carlo, cuya aplicacion se facilita cada dia mas debido al rapido avance de la tecnologia de 
computation. 

En este trabajo se desarrolla la estimation de un modelo estocastico de dicha variacion 
temporal de las ondas sismicas para la ciudad colombiana de Manizales, con base en los 
registros obtenidos alb en aiios recientes. Como aplicacion del modelo se presenta un analisis 
de Monte Carlo de una estructura representativa de la epoca en la que no se practicaba el 
diseno sismo-resistente en la ciudad, que corresponde a los anos anteriores a 1981. El objetivo 
es examinar la vulnerabilidad de dicho tipo de construction ante un sismo que tenga una 
probabilidad de excedencia de 10 % en 50 anos, lo que corresponde a un perfodo de retorno 
de, grosso inodo, 500 anos. 



2. REGISTROS SISMICOS OBTENIDOS EN MANIZALES 

La tabla 1 presenta los datos generates de los registros obtenidos en la ciudad que se 
ban utilizado para este estudio. Los cinco primeros fueron obtenidos en un acelerometro 
tipo Montana de registro fotografico en papcl, instalado en la ciudad por el U. S. Geological 
Survey y administrado por el Instituto Geofisico de los Andes. Cuatro de estos registros 
fueron digitados por Cordoba y Gomez (1987), y el quinto por Hurtado et al. (1981). Los 
registros restantes ban sido obtenidos en los equipos digitales marca Kinemetrics gestionados 
por el autor entre 1991 y 1993 en la Universidad Nacional en Manizales. 

Puede decirse que casi todos los registros corresponden al perfil tipico de los suelos 
de la ciudad, descrito por Arango y Gutierrez (1992), en el que los estratos superficiales 
estan formados por limos arenosos de origen volcanico y de gran espesor. A este perfil 
tipico pertenecen los sitios indicados en la tablas 1 como El Cable, Banco del Comercio, 
Universidad Nacional, Confamiliares y E. P. M. Las principales variaciones que pueden 
ballarse a lo largo y ancho del casco urbano corresponden a la presencia de algunos rcllenos 
mecanicos o bidraulicos de diferente espesor, los cuales alteran el comportamiento de las 
ondas sismicas en una medida que resulta dificil de ponderar. Las estaciones que figuran en 
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Tabla 1 Registros sismicos en Manizales 



Sismo No. 


Registros No. 


Epicentro 


Fecha 


Estacion 


i 


1 


Dabeiba 


08.31.77 


Banco del Cornercio 


2 


2 


Pereira 


04.13.75 


Banco del Cornercio 


3 


3 


Versalles 


06.25.80 


Banco del Cornercio 


4 


4 


Umpala 


03.22.77 


Banco del Cornercio 


5 


5 


La Tebaida 


05.18.76 


Banco del Cornercio 


6 


6, 7 


Belalcazar 


08.15.92 


El Cable 




8, 9 


Belalcazar 


08.15.92 


Universidad Nacional 


7 


10, 11 


Murindo 


10.17.92 


El Cable 




12, 13 


Murindo 


10.17.92 


Universidad Nacional 


8 


14, 15 


Murindo 


10.18.92 


El Cable 




16, 17 


Murindo 


10.18.92 


Universidad Nacional 


9 


18, 19 


Murindo 


10.18.92 


El Cable 




20, 21 


Murindo 


10.18.92 


Universidad Nacional 


10 


22, 23 


Toribio 


06.06.94 


J. Hada 




24, 25 


Toribi'o 


06.06.94 


E.P.M. 


11 


26, 27 


Tauramena 


01.19.95 


E.P.M. 


12 


28, 29 


Calirna 


02.08.95 


El Cable 




30, 31 


Calirna 


02.08.95 


Confamiliares 




32, 33 


Calima 


02.08.95 


E.P. M. 


13 


34, 35 


Risaralda 


08.19.95 


El Cable 




36, 37 


Risaralda 


08.19.95 


E. P. M. 




38, 39 


Risaralda 


08.19.95 


Acueducto 



las tablas como Universidad Nacional, El Cable y Confamiliares presentan un perfil similar, 
con la diferencia de que en el ultimo hay una capa de espesor moderado de dichos rellenos que, 
aparentemente, parece haber ocasionado pequenas amplificaciones en ciertos perfodos. For 
otra parte, el suelo del edificio E.P.M. y zonas aledanas presenta caracterlsticas excepcionales 
que ban causado tradicionalinente mayores intensidades alii que en otras partes de la ciudad. 
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Registro original 
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Figura 1 



Acelerograma del registro No. 7 (unidades: s, cm/s 2 ) 



3. ESTIMACION DE PARAMETROS DEL MODELO ESTOCASTICO 



La figura 1 muestra la historia de aceleraciones del regsitro No. 7. En ella pueden 
verse dos rasgos tipicos de la naturaleza no estacionaria de los acelerogramas, a saber, 
(a), las fuertes variaciones en amplitudes, que hacen que en general consistan de una zona 
ascendente, una fase de movimiento fuerte y una zona de desvanecimiento; y (b), la variacion 
del contenido frecuencial a lo largo del evento, asociado a las velocidades propias de cada tipo 
de ondas y a sus propias naturaleza / formacion. En consec’uencia, un modelo estocastico 
riguroso de la action sismica debe incluir su no estacionaridad en amplitud y en frecuencia. 
Para ello se ha adopt ado el modelo de espectro instantaneo (Yeh and Wen 1990), en el cual 
la densidad espectral de potencia evolutiva G e (o;,t) esta dada por 



(i) 

donde £(t) una funcion de variacion de las amplitudes y K,(t) una funcion que describe la 
evolution de las frecuencias dominantes del registro. Esta s funciones se estiman a partir de 
registros reales, lo que permite obtener los parametros que caracterizan el modelo global de 
un proceso estacionario, caracterizado por una densidad espectral G(uj). En este caso, el 
objetivo es identificar los parametros del modelo espectral de Clough-Penzien (Clough and 
Penzien 1993), dado por 
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GH = 



U) 



(«;-«*)*+ 4 «/> f v 



w ;+4,; w y 

(^~ u y + 

o 8 8 



(2) 



donde G w es la densidad espectral de un ruido bianco de fondo, u> g , i/ g son la frecuencia y el 
ainortiguamiento de un filtro que caracteriza basicamente al suelo regional (conocido como 
filtro de Kanai - Tajimi; cf. Kanai 1961; Tajimi 1960), mientras que w f y parametros 
de un filtro paso-alta. Los valores mas importantes del modelo son, con mucho, w , v , ya 
que Wpt'f solamente determinan la cafda de la densidad espectral hacia cero en la zona 
de muy bajas frecuencias (ver figura 2), mientras que G w puede determinarse a partir de 
la aceleracion maxima que se busque dar a los acelerogramas sinteticos. Notese que en la 
ecuacion (1) la division de u> por la derivada de la funcion de modulacion de frecuencias causa 
una modificacion temporal de la densidad espectral de la forma indicada esquematicamente 
por la figura 2. 




Figura 2 Densidad evolutiva tipo Clough-Penzien (unidades: 
rad/s, cm 2 /s 3 ) 



El proceso de calculo de los parametros que caracterizan del modelo es el siguiente: 

1. Calculo de las funciones empfricas de modulacion de amplitudes, £(t) y de frecuencias 
«(*). La primera se estima a partir de la energia del acelerograma, cuya ecuacion 
aparece mas adelante. La segunda se calcula a partir del valor acumulado de cruces 
del acelerograma por el nivel de aceleracion cero (cero- cruces). 

2. Estabilizacion estacionaria del registro, es decir, calculo de un registro estacionario equiv- 
alente obtemdo al remover las tendencias no estacionarias en amplitud y frecuencia car- 
acterizadas por £(t) y K (t). En este estudio se utilizaron las funciones empfricas de 
amplitud y frecuencia en lugar de modelos no lineales ajustados a ellas. Esto se debe 
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a que la ultima opcion tiene sentido cuando se busca generar acelerogramas sinteticos 
similares a uno dado (por ejemplo, en estudios de danos caudados por un evento que haya 
sido registrado) o cuando se pretende realizar un estudio estadistico de los parametros 
de las funciones £(t) y K,(t). Como en el caso presente no se trata de modelar un evento 
dado y, ademas, la base de datos disponible no es lo suficientemente amplia como para 
obtener una modelacion fiable de los multiples parametros de ambas funciones, se ha 
optado por utilizar las funciones empfricas directamente en la estabilizacion estacionaria 
de cada registro y concentrar el estudio estadistico en los parametros de la densidad 
espectral u; g y v % . 

Este proceso se realiza en dos fases. Inicialmente se estabiliza el registro en amplitudes 
dividendo la senal por la funcion de amplitudes empirica o ajustada: 



m, (t) 



a(t) 

m 



(3) 



Luego, la senal asi obtenida se mapea del eje de tiempo ficticio constituido por la funcion 
de modulacion de frecuencias al eje de tiempo real: 



m(t) = rrij («(<)) (4) 

Como ilustracion de este proceso, la figura 3 muestra la historia de aceleracion esta- 
bilizada correspondiente al registro No. 7 de la figura 1. Las densidades espectrales de 
potencia de los registros estabilizados se calcularon como el promedio de las densidades de 
tres segmentos de los rnismos. Para el calculo se utilizo la ventana espectral de Hanning 
(cf. Priestley 1981). 

3. Calculo de la duracion del segmento de fase fuerte, .s 0 , segun la definition de Vanmarcke 
y Lai (1980): 



y 



= >’ s ” i 1 - 367 '- 



5 =2 — 

0 max(a(^)) , 



s 0 < 1-36T. 



(5) 

( 6 ) 



donde T s es el periodo dominante de las ondas en la fase fuerte y E ^ es la energia del 
registro dad a por 



E a 




(7) 



Esta duracion se utilizara para modelar la modulacion de amplitudes en la generacion de 
acelerogramas sinteticos que se discute mas adelante. 



4. Calculo de los parametros espectrales del modelo de Clough y Penzien. A los parametros 
a ; f y se les asignaron unos valores fijos que se juzgaron adecuados a la vista de los 
espectros empiricos de los registros etabilizados, mientras que los parametros de la 
parte de Kanai - Tajimi del espectro se calcularon siguiendo el metodo de ajuste de 
los momentos espectrales (Lai 1982; Faravelli 1988a). Este metodo consiste en calcular 
los parametros del modelo a partir de la igualacion de los momentos espectrales de este a 
los obtenidos a partir del registro estacionario. Los momentos espectrales estan definidos 
como 
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Figura 3 Acelerograma estacionario obtenido del registro No. 7 
(unidades: s, cm/s 2 ) 



A j= uj J G(uj)duj 

J(\ 



( 8 ) 



Para calcular los tres valores desconocidos del modelo {G w , w g y i/ g ) se requiere resolver 
simultaneamente las ecuaciones correspondientes a los momentos de ordenes 0, 1 y 2 por 
medio de un algoritmo de estimation no lineal (Bard 1974; Press et al. 1992). Para ello 
se utilizaron las ecuaciones explicitas de los momentos espectrales del filtro de Kanai - 
Tajimi obtenidas por Faravelli (1988a). El algoritmo de ajuste no lineal utilizado fue el 
Levenverg - Marquart. Por otra parte, con el fin de evitar un sesgo en la estadfstica de 
los valores se prescindio de algunos registros del sismo No. 8 y de todos los del No. 9, 
que es una replica del anterior. 

La tabla 2 reune los valores de los parametros a 0 , w g y v % . En las figuras 4 y 5 se presentan 
histogramas de los dos ultimos valores junto con las funciones de densidad siguientes (tipos 
Weibull y Lognormal, respectivamente) que se han juzgado adecuadas para ellas: 



/K) = ^(^Ll^)^5 exp[ _ ( iiJL^M.55 1 



5.1 



5.1 



5.1 



-) 1 ' 50 ] 



( 9 ) 



Siy t ) 



1 r l.lni/ g + 1.9442 v2l 

0.3865^ g eXI>[ 2 0.3865 )] 



( 10 ) 
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Tabla 2 Valores calculados de los parametros del modelo estocastico 



Registro No. 


50 

(s) 


a max 

(cm/s 2 ) 


o; g 

(rad/s) 




5 


9.73 


48.60 


8.85 


0.150 


6 


12.04 


16.30 


16.34 


0.247 


7 


20.24 


16.25 


17.49 


0.262 


8 


14.31 


16.29 


13.93 


0.250 


9 


13.33 


14.88 


15.55 


0.233 


10 


21.34 


6.22 


10.92 


0.093 


11 


21.36 


6.75 


10.86 


0.135 


12 


22.17 


6.34 


10.17 


0.093 


13 


22.85 


5.86 


9.11 


0.119 


14 


30.70 


14.83 


12.09 


0.145 


15 


32.99 


13.92 


11.03 


0.149 


24 


21.46 


9.44 


8.83 


0.052 


25 


29.93 


5.90 


9.24 


0.069 


26 


20.34 


12.30 


9.11 


0.077 


27 


13.16 


19.99 


9.32 


0.086 


28 


16.53 


26.42 


13.66 


0.191 


29 


20.46 


38.02 


13.48 


0.184 


30 


25.34 


22.98 


11.23 


0.133 


31 


14.10 


28.53 


11.31 


0.122 


32 


12.94 


55.86 


9.01 


0.136 


33 


25.46 


35.38 


9.53 


0.110 


34 


16.14 


51.58 


17.03 


0.229 


35 


11.91 


46.78 


18.93 


0.231 


36 


19.30 


35.26 


13.22 


0.181 


37 


24.46 


34.06 


12.16 


0.178 



La figura 6 muestra la densidad espectral evolutiva del registro 7 correspondiente al 
modelo de Yeh-Wen, la cual fue calculada con los parametros que aparecen en la tabla 
2 y usando sus funcion empirica de evolucion de frecuencias k(t). Como modulacidn de 
amplitudes se uso la funcion de Amin y Ang (1966) 

* < *, 

= < i, t,<t<t 2 

{ t 2 < t 

con parametros t x = 2 s, c = 0.18 y t 2 = 4- s 0 . 



(ii) 
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Figura 4 Histograma y modelo probabilista de la frecuencia de 
Kauai- Tajimi para Manizales (unidades: rad/s) 




v 

9 



Figura 5 Histograma y modelo probabilista del amort iguamiento 
de.Kanai- Tajimi para Manizales 
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Figura 6 Densidad espectral evolutiva basada en el registro No. 7 
(unidades: s, cm/s 2 ) 



4. GENERACION DE ACELEROGRAMAS ARTIFICIALES 

Con el fin de utilizar el modelo espectral asi definido para generar acelerogramas arti- 
ficiales adecuados para la ciudad consistentes con su amenaza sismica, se hace necesario 
examinar la relacion existente entre los valores s 0 , u; g y v de un lado y la aceleracion 
maxima, a max , por otro, debido a que la amenaza sismica se encuentra definida para todo 
el territorio nacional de Colombia enlerminos de esta (Garcia et al . , 1984). La figura 7 
muestra la relacion existente entre s 0 y a max . En ella tambien se recogen los datos de la 
costa oeste norteamericana evaluados por Lai (1982). Puede observarse que los datos de 
Manizales muestran la misma tendencia a una correlation negativa de s 0 con respecto a 
flmax) lo cual se explica principalmente por el hecho de que la duracion crece con la distancia 
epicentral al contrario que la aceleracion maxima. Asimismo, se puede ver que los datos 
de Manizales, en el pequeiio rango de aceleraciones que ha sido posible registrar hasta el 
momento, muestran una dispersion muy inferior que los usados por Lai (1982), lo que se 
debe a la mayor homogeneidad de la muestra de Manizales. Esto ilustra la inconveniencia 
de utilizar informacion de bases de datos correspondientes a regiones diversas en estudios de 
objetivo local como aquellos para los cuales un trabajo como este pretende servir de insumo, 
es decir, la generacion de acelerogramas artificales para estudios de vulnerabilidad y riesgo. 

Sobre la base de la informacion reunida en la tabla 2 se ha calculado la siguiente regresion 
entre duracion de la fase fuerte y aceleracion maxima: 
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Figura 7 Relacion aceleracion maxima / duracion de la fase fuerte 
(unidades: s, cm/s 2 ) 



s 0 = -0.1849 a max + 24.0612 + € (12) 

En csta expresion e es una variable aleatoria normal con media cero y desviacion estandar 
5.58 s. Como quiera que para valores altos de aceleracion se puede obtener duraciones 
negativas, se propone fijar un valor mfnimo de 1 s. La relacion se postula lineal pues la 
escasez de datos no permite suponer otra cosa. 

De otra parte, las figuras 8 y 9 muestran las relaciones u> g - a max y v - a max , respec- 
tivamente. En consonancia con lo oLservado por Lai (1982), ninguna de las dos variables 
muestra una tendencia definida a crecer o disminuir con la aceleracion maxima, lo cual hace 
que, para efcctos practicos, puedan ser consideradas efectivamente como propiedades (es- 
tocasticas) de los suelos de la region bajo estudio, modeladas por las ecuaciones (8) y (9) 
rcspectivamente. Para estudios de vulnerabilidad y riesgo, en los cuales se hace necesario 
generar acelerogramas artificiales para un rango amplio de aceleraciones maximas, esta in- 
dependencia de los parametros espectrales del nivel de aceleracion sismica permite generar 
los acelerogramas a partir de una misma densidad espectral caracterizada ella misma por 
parametros aleatorios. En consecuencia, el proceso de generacion de acelerogramas artifi- 
ciales propuesto para la ciudad, a partir del modelo estocastico asf definido, es el siguiente: 

1 . Definir la aceleracion maxima objetivo del registro. 

2. Generar una duracion aleatoria de la fase fuerte, teniendo en cuenta que el parametro e 
en la ecuacion (12) es aleatorio. 
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3. Generar valores aleatorios de los parametros de la densidad espectral de potencia de 
acuerdo a sus distribuciones (ecuaciones 9 y 10). 

4. Generation de un ruido bianco de duracion equivalente a la duracion total del evento y 
aplicacion de la funcion de modulacion de Amin y Ang (ecuacion 11). 

5. Simulation del acelerograma por medio de la solution de para U [ de la dinamica del filtro 
variable de Clough - Penzien: 

U % + {2v % u g k - = -k 2 at)W(*(t)) (13) 

K> 

U, + (2 u f u) f k - %)U f + uj*kU, = -2 v g u g kU g - u*KU g (14) 



6. Correction de la lfnea de base. 



5. VULNERABILIDAD DE EDIFICIOS SIN DISENO SISMO-RESISTENTE 



Como una primera aplicacion practica del modelo estocastico expuesto, en este apartado 
se hara uso de el para estimar la vulnerabilidad de una estructura representativa de las 
const rucciones de concreto reforzado carentes de diseno sismo-resistente bajo condiciones 
de un sismo fuerte. Como valor indicativo de este ultimo se toma el correspondiente a una 
probabilidad de excedencia de 10% en 50 ahos, caracterizado por una celeracion maxima 
de 0.25 g y reglamentado por el codigo de diseno sismo-resistente de la ciudad para las 
construcciones nuevas. La meta de este analisis es la obtencion de una funcion de distiibucion 
de probabilidad del daho de una estructura representativa de los edificios de la ciudad de la 
epoca anterior a la expedition de la primera norma sismica (1981), la cual permite estimar 
el escenario posible de daiios y perdidas en ese grupo de edificios. 

El metodo que se seguira para esta evaluation es el metodo de Monte Carlo (Hurtado and 
Barbat 1998), que consiste en la realization de un arnplio numero de analisis deterministas 
de la estructura con datos aleatorios con el fin de disponer de una muestra sobre la cual 
realizar las estadisticas correspondientes. Los pasos requeridos por el metodo en un caso 
como este son los siguientes: 

1. Diseno de una estructura altamente representativa de hi epoca. La estructura en cuestion 
es un edificio de concreto reforzado de seis plantas, con cinco marcos de dos vanos cada 
uno, que soportan losas de concreto y cargas vivas de tipo vivienda. El diseno se realizo 
de acuerdo a las normas y costumbres de las decadas de los sesenta y setenta. 

2. Definicion de un indice de dano. En los liltimos anos se ha impuesto en la literatura 
internacional el modelo de Park y Ang (Park 1984; Ang 1987) el cual estima el dano 
producido en cada elemento estructural como funcion de la deformation maxima y de la 
energia disipada. La ecuacion del indice para un elemento i es 



d '~ S f + wJ iE 



(15) 



donde y son la deformation maxima bajo carga sismica y bajo carga monotonica. 
respectivamente, F y la resistencia a cedencia del elemento y (3 un coeficiente de peso de 
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Figura 8 Relation aceleracion maxima/ frecuencia del filtro (uni- 
dades: cm/s 2 , rad/s) 
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Figura 9 Relacion aceleracion maxima/ amortiguamiento del filtro 
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la energia disipada por histeresis, E. Park (1984) propone expresiones para las variables 
mencionadas y modelos de histeresis para cada tipo de carga a los que tipicamente se ven 
sometidos los elementos estructurales bajo sismos: flexo-compresion, corte y adherencia. 
Asimismo, para un edificio en su conjunto, propone tomar como indice de daho global la 
expresion 



Z Eidi 






(16) 



que implica el tomar la energia disipada como factor de ponderacion. La mayor ventaja 
de este indice es que ha sido cuidadosamente calibrado con respecto a dahos sismicos 
realmente ocurridos por medio de analisis no lineales de historia de respuesta de varios 
edificios sometidos a los registros de los terrcmotos correspondientes (Ang 1987). La 
tabla 3 muestra los rangos de valores del indice correspondientes a varias calificaciones 
del daiio. los cuales suelen ser tornados como referenda en diferentes estudios de este tipo 
(Barbat et al. 1995; Singhal and Kiremidjian 1996; Abrams and Shinozuka 1997). 



Tabla 3 Rangos del indice de Park y Ang para diferentes estados de 
daiio (Singhal and Kiremidjian 1996) 



Estado de daho 


Rango 


Leve 


0.1 - 


0.2 


Moderado 


0.2 - 


0.5 


Fuerte 


0.5 - 


1.0 


Colapso 


> 1 


.0 



3. Definicion de las variables aleatorias. En este caso se tomaran las variables que se indican 
en la tabla 4. La variable /' es la resistencia ultima del concreto a compresion, cuya media 
se toma igual a 1.14 veces la resistencia nominal del concreto de diseiio cornun en la epoca, 
0.02058 kN/mm 2 . En cuanto a la media de la resistencia ultima del acero / , se toma un 
cinco por ciento superior a la nominal, 0.4116 kN/mm 2 . 

4. Generacion de conjuntos de datos aleatorios y combination de los mismos de acuerdo a 
una tecnica de reduction de la varianza, que para este caso sera la de muestreo descriptivo 
(Ziha 1995). 

5. Generacion de un conjunto igual de acelerogramas sinteticos por medio de la solution de 
las ecuaciones del filtro variable de Clough-Penzien 

6. Analisis de historia de respuesta no lineal de todos los modelos con los grupos de datos 
aleatorios. Para este proposito se utilizo el programa IDARC (Kurmath et al. 1992), 
desarrollado para el analisis inelastico de estructuras de concreto reforzado. 

7. Analisis estadistico de los resultados. 

Se realizaron en total 66 analisis inelasticos del modelo. La figura 10 muestra los valores 
de la distribution de probabilidad empirica del mdice de daho. A ellos se ha ajustado 
una funcion lognormal, Ja cual permite estimar las probabilidades asociadas a estados 
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Tabla 4 Descripcion de las variables aleatorias 



No. 


Variable 


Distribucion 




a 


1 


e 


Normal 


0 


5.58 


2 


Ss 


Weibull 


12.096 


3.022 


3 




Lognormal 


0.154 


0.062 


4 


/ 


Normal 


0.0239 


0.0033.5 


5 


4 


Lognormal 


0.4410 


0.04851 



Tabla 5 Probabilidades de excedencia del mdice de dano 



d 


p[d >d\ 


0.1 


9.998 E -01 


0.2 


5.060 E -02 


0.3 


1.409 E -06 


0.4 


3.864 E -12 




Figura 10 Funciones de distribucion del mdice de dano de Park v 
Ang 
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no observados en la simulation, indicados en la tabla 5. Alii puede observarse que para 
el grado de daiio leve hay casi total certeza de su ocurrencia en el tipo de estructuras 
consideradas, mientras que las probabilidades correspondientes a niveles de daiio mayores 
son muy reducidas. El fuerte contraste entre estos valores de probabilidad de excedencia 
entre los grados leve y moderado se explica por el reducido valor del coeficiente de variation 
del indice de daiio, que es tan solo de 0.133. Una dispersion menor aun se obtuvo en el caso 
del radio de deriva maxima (definida como la relation entre la deriva de piso y su altura). 
La media y el coeficiente de variacion de esta variable fueron, respectivamente, 0.0092 y 
0.05. Con el fin de valorar este resultado, es conveniente mencionar que los codigos de 
diseiio sismo-resistente especifican un valor maximo del radio de deriva entre 0.01 y 0.015 
como medida de prevention del colapso. Esto, junto a la valoracion probabilista del indice 
de daiio, permite concluir que, al menos teoricamente, las estructuras regulares sin diseho 
sismico de la ciudad no presentan una vulnerabilidad alta ante los sismos mas fuertes que 
se esperan en ella y que es poco probable que sufran dahos considerables en tal caso. Con 
el fin de complementar esta diagnostico preliminar deben efectuarse estudios similares en 
estructuras irregulares sin diseho sismo-resistente. 



CONCLUSIONES 

En este trabajo se ha expuesto la estimation de un modelo estocastico de la action 
sismica sobre una ciudad colombiana, que permite realizar diferentes estudios probabilistas 
de vulnerabilidad y riesgo de estructuras por medio de metodos analiticos o por tecnicas de 
Monte Carlo. De las multiples variables que componen el modelo evolutivo se ha centrado 
la atencion en la duration de la fase fuerte y en la frecuencia y amort iguamiento del filtro, 
los cuales se han estimado por tecnicas de parametrizacion no lineal. Las funciones de 
probabilidad y las regresiones propuestas sirven como base para generar los acelerogramas 
artificiales requeridos para los estudios de vulnerabilidad sismica por el metodo de Monte 
Carlo. Una aplicacion desarrollada en este sentido muestra que los edificios regulares 
construidos en la ciudad sin diseho sismo-resistente presentan una vulnerabilidad reducida 
ante movimientos considerados como fuertes en las normas de construction de la ciudad. 
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RESUMEN 

En el presente artfculo se describe un sistema de refuerzo de estructuras de hormigon aplicando 
materiales compuestos, constituidos por fibras de carbono y resina epoxidica, en forma de 
laminados adheridos extemamente a la estructura. Junto a una breve description del concepto de 
materiales compuestos, propiedades y aplicaciones de los mismos, se expondran las principales 
razones que han convertido a este sistema de refuerzo en una autentica altemativa frente a metodos 
convencionales, principalmente a flexion. Posted ormente se comentaran los criterios de diseno 
empleados en el proyecto y diseno de las estructuras reforzadas que se presentan (particularizados, 
fundamentalmente, para los laminados unidireccionales constituidos por fibras de carbono y resina 
epoxidica), asi como los diversos procesos constructivos que se llevan a cabo segun el tipo de 
material compuesto aplicado en el refuerzo. A1 final del trabajo los autores exponen las principales 
conclusiones que han obtenido a traves de los refuerzos con materiales compuestos con fibra de 
carbono realizados. 



1 INTRODUCTION 



La aplicacion de laminados compuestos de fibra de carbono y resina epoxidica en el refuerzo de 
estructuras de hormigon, metalicas, de madera, de piedra, etc. comienza a constituir, en la presente 
decada, una autentica altemativa al sistema de refuerzo convencional median te encolado de chapas 
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de acero (fundamental rnente en refuerzos a flexion), tanto por sus mejores prestaciones mecanicas 
especificas y resistencia a la corrosion, como por los ahorros obtenidos en el proceso total del 
refuerzo, gracias a la ligereza que presentan los materiales compuestos, facilitando su transporte, 
manejo y puesta en obra, empleando medios auxiliares ligeros durante cortos periodos de tiempo. 

La primera aplicacion en Espana del presente sistema de refuerzo se llevo a cabo en el puente del 
Drago (Barcelona, 1.996). A nivel mundial, la primera estructura mediante laminas encoladas de 
materiales compuestos con fibra de carbono data de 1 .99 1 : el puente Ibach, en Lucema. 

Tan solo en Suiza y Aleman ia, el numero de realizaciones de refuerzos de estructuras (tanto 
puentes como edificios) con laminas compuestas encoladas asciende a mas de doscienta cincuenta, 
existiendo mas de un miliar mundialmente (concentradas, fundamental mente, en Suiza, Alemania, 
Japon, EE.UU. y Canada). En Espana, en la actualidad, podemos encontrar varias decenas de 
realizaciones (1). 

2 MATERIALES COMPUESTOS 

El termino material compuesto se reserva para aquellos materiales bifasicos (fibra + matriz) 
fabricados expresamente para mejorar los valores de las propiedades que los materiales 
constituyentes presentan por separado, exibiendo una interfaz identificable entre ellos. 

Si bien la variedad de los mismos es enorme, se cumple como denominador comun la existencia 
de un elemento fibroso, que aporta rigidez y resistencia, con un ratio longitud/espesor claro, y una 
matriz (generalmente resina) que configura geometricamente el material compuesto. 

Normalmente, la matriz es flexible y poco resistente, con la mision principal de transmitir los 
esfuerzos de unas fibras a otras, y entre ellas y la superficie adyacente, ademas de proteger a las 
fibras de posibles danos mecanicos y ambientales. 

Las fibras empleadas pueden estar constituidas por carbono, vidrio, boro, aramida, metal, 
ceramica, etc. La naturaleza de las matrices tambien puede ser muy variada: organicas, minerales, 
metalicas, ceramicas, etc. Referente a las aplicaciones en la ingenieria civil, el mercado de los 
materiales compuestos lo conforman, mayoritariamente, las matrices organicas y las fibras de 
carbono, vidrio y aramida. 

2.1 Propiedades de los materiales compuestos 

Como propiedades generales podemos destacar: elevadas resistencias mecanicas 

(comportamiento lineal hasta rotura), resistencia a la corrosion, resistencia a los ataque quimicos, 
ligereza, presentacion en cualquier longitud (eliminacion de juntas), muy buen comportamiento a 
fatiga, buenas propiedades dielectricas, etc. 

Cuando las fibras, dentro del material compuesto, se encuentran orientadas en una unica 
direccion (presentacion industrial conocida como roving ), nos encontramos con un material 
anisotropo, obteniendose la maxima resistencia y modulo en la direccion del eje de las fibras. Una 
disposicion del refuerzo en dos direcciones (presentacion denominada tejido) produce diferentes 
resistencias segun el angulo que formen las fibras. En el caso de que las fibras se encuentren 
orientadas en todas direcciones (presentacion denominada mats), el comportamiento del compuesto 
es isotropico, reduciendose sustancialmente los valores de la resistencia en relacion con los 
obtenidos en los compuestos tipo roving segun la direccion de las fibras. 

La forma de los materiales compuestos mas usual empleada en aplicaciones estructurales se 
denomina laminado, constituido por varias laminas de fino espesor. Dependiendo de la orientacion 
de las fibras en cada lamina y de la secuencia de las mismas dentro del laminado se puede generar 
un amplio rango de propiedades mecanicas y ffsicas. 

Las propiedades mecanicas- del material compuesto, en cualquier direccion, dependen del 
porcentaje de fibras, en volumen, orientadas en la misma direccion, asi como de las propiedades 
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mecamcas de las fibras, de la longitud. forma y composicion de las fibras, de las propiedades 
mecamcas de la resina y de la adherencia entre las fibras y matriz. Dichas propiedades dependen 
sobremanera de la direccion de medida en relacion con la direccion de las fibras. 

Las propiedades de la matriz influyen en la resistencia a cortante del material compuesto, asi 
como en las propiedades quimicas, electricas y termicas del compuesto, ademas de dar soporte 
lateral contra el pandeo de las fibras bajo solicitaciones de compresion (2). 

Para poder apreciar esta variabilidad de las propiedades de los materiales compuestos en funcion 
del tipo de ftbra y de resina empleados, se adjunta en la figura n° 1 la relacion tension - deformacion 
(ensayo uniaxial) para diversas fibras, asi como para el acero de armado y pretensado; 
constatandose el caracter completamente lineal entre tensiones - deformaciones de los materiales 
compuestos llevados hasta rotura. 
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* igura 1 Densidadcs y curvas tension - dcformacidn (ensayo uniaxial), para 
diversas fibras sintcticas y aceros de armado y pretensado (5). 

Dentro de las caracterfsticas mas importantes que presentan algunos materiales compuestos 
liente a otros tenemos: excelente adherencia fibra - matriz (vidrio y carbono), estabilidad 
dimensional (bajo coeficiente de dilatacion) (vidrio y carbono), incombustibilidad (vidrio). 
imputrescibilidad (vidrio), buenas propiedades dielectricas (vidrio y carbono), etc. 



2.2 Aplicaciones de los materiales compuestos 



En todos aquellos casos donde la ligereza, la resistencia mecanica especifica, la resistencia a la 
coirosion, el aislamiento electrico, la estabilidad dimensional, etc. sean importantes, se puede hablar 
de una presente o potencial apiicacion de los materiales compuestos. 

Dichas aplicaciones las encontramos en la industria automovilistica, naval, aeronautica, 
eportiva, electrica, defensa, construccion, bienes de consumo, materiales industriales y agricolas 
etc. Sorprende que, a pesar del impacto que tienen los materiales compuestos en otras areas, su 
presencia en la Ingenieria Civil ha sido virtualmente inexistente. 

No obstante, las aplicaciones llevadas a cabo en Construccion las encontramos en: estructuras de 
horrrugon en amb.entes agresivos, plataformas Off - Shore, depositos (la primera apiicacion 
estructural de los modemos materiales compuestos fue la construccion de depositos a presion 
mediante enrrollamiento continuo de fibras de vidrio), anclajes al terreno, construcciones no 
conductivas, no magneticas (radar, etc.), tuberias, refuerzos de estructuras, armaduras pasivas, 
armaduras activas, cables, paneles (edificacion), tableros de pasarelas (pasarela de Aberfeldy, Reino 
Unido), perfileria, recubrimientos de tuneles, etc. 

Debido a la enorme potencialidad, aplicabilidad y variedad de los materiales compuestos se 
requiere un detallado y juicioso estudio de seleccion de los mismos con objeto de conseguir una 
relacion cal idad/coste optima y efecti vidad estructural . 
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3 REFUERZOS DE ESTRUCTURAS DE HORMIGON 

Diferentes criterios pueden llevamos a refozar estructuras existentes: criterio de rigidez 
(reduccion de flechas en servicio), criterio de resistencia (mejora en la capacidad maxima de carga) 
y criterio de durabilidad (limitaciones en el ancho y distribution de fisuras en el hormigon). 

Las causas que pueden provocar dichas necesidades de refuerzo son muy diversas: dafios 
mecanicos en las estructuras provocados por accidentes o impactos de vehiculos, cambios 
funcionales, errores de calculo en el proyecto o durante la ejecucion, incremento de las cargas de 
trafico, corrosion del refuerzo (acero), etc. 

De entre los diversos metodos empleados en el refuerzo de estructuras, el encolado de chapas de 
acero mediante adhesivo estructural constituye uno de los sistemas mas ampliamente utilizado y de 
mayor eficacia (fundamentalmente en zonas traccionadas de refuerzos a flexion), sobre todo desde 
el desarrollo experimentado en los ultimos anos en las tecnicas adhesivas (la union de la chapa de 
acero al sustrato mediante adhesivo estructural ha demostrado ser mas eficaz que las realizadas con 
roblonado, atomillado, etc., gracias a una transmision mas uniforme de esfuerzos obtenida a traves 
de la interfaz laminado - sustrato), (el primer caso de aplicacion en edification se produjo en 
Durban, Sudafrica, en 1 .964 y el primero en puentes se llevo a cabo en Swanley y Kent - Inglaterra- 
en 1.975-77) (6). 

A pesar de la eficacia demostrada en el refuerzo mediante laminas encoladas de acero, nos 
encontramos con los siguientes incovenientes manifiestos: 

• Las laminas de acero constituyen elementos pesados en su manejo, transporte y 
ejecucion, limitando las longitudes maximas de las chapas a 6 - 8 metros (surgiendo la 
problematica de las juntas). Debido a su elevada densidad, requieren medios auxiliares 
importantes para su apuntalamiento, empleados durante largos plazos de tiempo, 
pudiendo provocar molestias por los cortes de trafico producidos o por la ejecucion de 
obras en edificios en funcionamiento. 

• Corrosion de la lamina de acero (afectando a la adherencia entre hormigon y chapa). 

• Debido a la gran rigidez que presentan las chapas de acero, se requiere una elevada 
planeidad en las superficies de encolado. 

Con la necesidad de solventar las deficiencias presentadas en la utilization de chapas metalicas, 
y ante el elevado y creciente numero de estructuras existentes que requieren ser reforzadas, se 
desarrollo en los Laboratories Federales Suizos para el Ensayo de Materiales (EMPA), a mediados 
de los anos ochenta, un sistema alternative de refuerzo a flexion a base de materiales compuestos 
(mas del 5% de los puentes deteriorados en Europa pueden ser reforzados mediante esta nueva 
tecnologia). En la actualidad, existen diversos laminados compuestos en el mercado para el refuerzo 
de estructuras. 

La resistencia a la corrosion de los laminados compuestos permite que la adherencia entre 
sustrato (hormigon, madera, etc.) y lamina no sufra alteraciones, asi como su reducida densidad 
frente al acero permite una mayor maniobrabilidad y facilidad en la ejecucion, requiriendose medios 
auxiliares ligeros con plazos de puesta en obra muy reducidos (las densidades de los diversos 
materiales compuestos pueden apreciarse en la figura 1). Esta mayor ligereza de los materiales 
compuestos permite su transporte en rollos de hasta 300 m de longitud, evitando la problematica de 
las juntas (figura 2). Como inconvenientes mas representatives podemos enunciar el alto precio del 
material y los posibles modos de rotura fragiles de no proyectarse adecuadamente el refuerzo. 
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Figura 2 Vista de la laminas de material 
compuesto (roving), enrolladas. 



Gracias al empleo de medios auxiliares mas ligeros y durante un plazo de tiempo menor, se 
pueden ilegar a obtener ahorros de hasta un 25% en el proceso total del refuerzo, compensando, de 
este modo, el mayor precio de los materiales compuestos frente al del acero. 

Con lo expresado anteriormente, se puede afirmar que esta reciente tecnologia de refuerzo a 
flexion mediante la aplicacion de materiales compuestos tiene un rango de aplicacion muy amplio, 
pudiendose llevar a cabo no solo donde los requerimientos de altas prestaciones, la maniobrabilidad 
y rapidez de la ejecucion constituyan un punto importante, sino tambien donde el planteamiento sea 
estrictamente economico, gracias al ahorro que se obtiene por la utilization de medios auxiliares 
ligeros con plazos de ejecucion rnmimos. 



4 CRITERIOS DE DISENO 



Como recomendaciones generales en el diseno de refuerzos de estructuras, tanto con materiales 
compuestos como con materiales mas conventional, podemos destacar las siguientes: 

* 

• Evitar el colapso de la estructura en el caso de elimination del refuerzo (accidente, 
vandalismo, etc.). 

• Contemplar el estado previo de deformaciones en la estructura antes de reforzar 
(redistribution de esfuerzos a nivel de estructura y seccional). 

• Estudio de la idoneidad del refuerzo seleccionado. 

El esquema general de calculo comprendena la comprobacion de la seguridad de la estructura sin 
reforzar (flexion, cortante, etc.), verificacion de la seguridad de la estructura reforzada y 
verificacion de los Estados Limites de Servicio. 

La fiabilidad de una estructura reforzada con materiales compuestos depende sobremanera del 
tip° de material compuesto empleado (las caractensticas meanicas de los mismos pueden variar con 
el tiempo, debido a fenomenos tales como el envejecimiento, etc.), la solicitation de trabajo, 
condiciones ambientales a las que se vera sometida, etc., resultando compleja la extrapolacion de los 
criterios de diseno y de los coeficientes de seguridad para los diferentes materiales compuestos, 
requiriendose un estudio particular y riguroso para cada tipo de material y aplicacion del mismo. 



80 



M.D.G. Pulido & J. A. Sobrino 



En la figura 3 se pueden apreciar los di versos modos de rotura que se pueden experimentar en 
elementos estructurales reforzados con materiales compuestos adheridos en zonas traccionadas 
(refuerzos a flexion). 




1: Rotura a traccion del laminado. 

2: Rotura por compresion en el hormigon. 

3: Plastificacion del acero de armado. 

4: Rotura a traccion del hormigon. 

5: Rotura interlaminar del laminado. 

6: Rotura de la pelicula adhesiva. 

7: Rotura adherencia (superftcie CFRP). 

8: Rotura adherencia (superftcie hormigon). 



Figura 3 Seccidn longitudinal de una viga reforzada con l£minas 
de materiales compuestos. Modos de rotura (6). 



Los danos aquf senalados pueden conducir a roturas inmediatas del elemento o a roturas locales 
con efectos negativos para su capacidad portante. 

En el Centro de Investigation suizo EMPA se han llevado a cabo (desde 1.985) de modo 
sistematico ensayos de estructuras reforzadas con laminas de materiales compuestos con fibra de 
carbono, bajo solicitaciones estaticas y dinamicas y con diferentes condiciones climaticas (7). Los 
modos de rotura observados fueron los siguientes: 

• Rotura de la lamina (tipo 1 en la figura 3) al alcanzarse la resistencia a traccion de la 

misma. " 

• Rotura en la zona de compresion del hormigon al superarse la deformacion maxima a 
compresion (tipo 2). 

• Rotura del acero al superarse su resistencia a fatiga o al excederse la deformacion de 
fluencia (tipo 3). 

• Desprendimiento de la lamina de material compuesto (tipo 4). Debido a un desplazamiento 
vertical de la superhcie inferior de la viga provocado por irregularidades locales por 
insuficiente preparation de la superficie y/o a desplazamientos relativos de las paredes de 
las fisuras por fuerzas cortante, se puede producir la deflexion de la lamina; provocando 
grandes esfuerzos de traccion perpendicular a la lamina que podrfan fisurar horizontal mente 
el hormigon y desprender la lamina en toda su longitud. 

• No se presento ninguna rotura del adhesivo a temperatura ambiente. 

Las roturas tipo 1 , 2 y 4 son graves por la rotura subita que se puede producir en el elemento, 
mientras que los modos 5 a 8 pueden provocar una rotura local que, eventual mente, puede ser 
reparada. 
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4.1 Dimensionamiento a flexion 



Hn el diseno de refuerzos a flexion se consideran, fundamentalmente, los mismos principios que 
se emplean en refuerzos mediante laminas de acero: 

• Equilibrio de la seccion. 

• Compatibilidad de deformaciones. 

• Hipotesis Navier - Bemouilli (secciones planas permanecen planas despues de la 
deformacion). 

N 4 ^°, S m<x * elos constitutivos del hormigon y el acero son los comunmente empleados (Codigo 
odelo, Eurocodigos, etc.), mientras que para el material compuesto se utiliza una relacion 
totalmente lineal entre tensiones y deformaciones hasta rotura. 

Si bien se trata de un metodo simplificado, los resultados obtenidos concuerdan con suficiente 
exactitud con los resultados experimentales. 

Con objeto de obtener una respuesta seccional ductil, asf como una mayor resistencia a flexion 
se debe d.senar el refuerzo de modo que la rotura de la lamina de material compuesto se produzca 
durante la plastificacion del acero y antes de que el hormigon alcance la rotura por compresion. En 
la hgura 4 se puede apreciar el comportamiento ductil a flexion de la seccion critica de la viga 
reforzada en el puente del Drago (diagrama momento - curvatura). 




Figura 4 Diagrama momento - curvatura, en $ecci6n critica de la viga 
retorzada con laminas de materiales compuestos con fibra de 
carbono. 



Para el caso concrete de refuerzos de estructuras mediante laminados unidireccionales (tipo 
rovings) de fibra de carbono con resina epoxidica, se exponen a continuacidn los principles 
entenos empleados en el diseno del refuerzo, basados en experimentos realizados en los 
laboratonos del EMPA (7), (8), (9), (10), (1 1), (12) y (13): 



Los anchos efectivos de la seccion, donde se puede considerar la colaboracion del 
laminado compuesto con el resto de la estructura, se detallan en la figura 5 



82 



M.D.G. Pulido & J. A. Sobrino 





MrO Mr. 



Mr= Mro+Mri 



Figura 5 Determinacion de anchos eficaces en secciones reforzadas. 



donde 

h: canto de la estructura. 

x: bloque de compresiones. 

s L : separacion entre ejes de laminados. 

b L : ancho del laminado. 

M r : Momento resistente de una seccidn de ancho s L 

M ro : Momento resistente de la zona reforzada. 

M R1 : Momento resistente de la zona no reforzada. 

• Con objeto de considerar el efecto de la fisuracion del hormigon se define un coeficiente de 
union, K, que relaciona los valores medios de las deformaciones con los valores maximos en el 
acero y laminado, afectando a la respuesta ultima seccional. Como se detalla en la figura 6, 
mientras en el hormigon fisurado a traccion se puede considerar constante la deformacion en el 
hormigon, las deformaciones en el acero y en los laminados compuestos oscilan alrededor de un 
valor medio (zona entre fisuras), exp^rimentando los valores maximos en las zonas proximas a 
las fisuras, traduciendose en puntas en las tensiones de trabajo. Como los diagramas momento - 
curvatura representan la respuesta estructural de la seccion en valores medios, la estructura no 
se agotara cuando se alcancen dichos valores, sino antes, debido al incremento de tensiones que 
se produce en las fisuras (el piano de deformaciones muestra los valores medios de las 
deformaciones en el acero y laminado, mientras que las tensiones y esfuerzos se refieren a la 
seccion fisurada, donde se producen los valores maximos). Los diferentes valores de k se 
expresan a continuation: 
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Figura 6 Distribution de deformaciones en zona fisurada. 



4.2 Dimensionamiento por deslaminacion 



Como se describio en el modo de rotura tipo 4 de la figura 3, el colapso de la estructura puede 
venir determinado por el desprendimiento del laminado (deslaminacion). En la figura 7 se detalla 
graficamente el proceso de deslaminacion: las fisuras (de ancho w) producidas por solicitaciones de 
flexion, en zonas traccionadas, pueden experimentar un desplazamiento vertical relativo, 8, entre los 
labios de la fisura por solicitaciones de cortante, sufriendo el laminado ese mismo desplazamiento 
vertical que, al estar traccionado, experimenta una fuerza de traccion perpendicular a la lamina, 
pudiendo provocar en el hormigon una fisura horizontal que se puede extender hasta el comienzo de 
la misma, desprendiendola. 




a) Aparicion de fisura 
por flexion. 




b) Desplazamiento vertical 
relativo entrc labios de 
fisura por cortante. 



c) Deslaminacidn del 
hormigon al superarse 
el cortante critico. 



Figura 7 Proceso de deslaminacion. 
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A continuation se adjunta una formula empiTica propuesta para determinar la fuerza cortante que 
puede provocar el desprendimiento de la lamina (fenomeno de deslaminacion), obtenida a partir de 
ensayos realizados en los laboratorios suizos EMPA, haciendo notar que, como en el caso del 
apartado de dimensionamiento a flexion, dichos resultados solo son aplicables para el caso concreto 
de laminados unidireccionales con fibra de carbono y resina epoxfdica, con un 60% en volumen de 
libra. Dicha formula solo es aplicable en elementos sometidos a carga uniforme, no siendo valida en 
casos donde se produzcan simutaneamente grandes solicitaciones de momentos flectores y esfuerzos 
cortantes, tales como elementos con grandes cargas puntuales o zonas de momentos negativos con 
grandes areas de cortante (voladizos y zonas alrededor de column&s). 



V A= T cA bx + *cA b ( h - x )( l -£u,,'Z)+ x(E s A s + E l A l ) con (l-£ 7j „/£)>0 



donde 

V A : cortante critico de deslaminacion. 

t ca : Tension cortante ultima del hormigon. 

b: Ancho de la viga o elemento. 

h: Canto de la viga o elemento. 

x: Posicion de la fibra neutra. 

Elov Deformacion de la lamina. 

E S A S : Modulo de deformacion y area del acero de refuerzo. 

ElAl* Modulo de deformacion y area del laminado. 

%, Constantes determinadas experimentalmente donde se tiene 
en cuenta la calidad del hormigon. 

Los valores de x y El™ se deben calcular en la section con esfuerzo cortante maximo. 

El primer sumando de la ecuacion que determina el valor de V A representa la transferencia de 
cortante en la zona de compresion del hormigon, el segundo la transferencia de cortante por la 
interfaz entre las paredes de la fisura (efecto de encaje) y el tercero la transferencia de cortante 
debido a las fuerzas verticales en el acero de refuerzo y en el laminado provocadas por el 
desplazamiento vertical 8 (efecto de armadura pasante). Todo ello se puede apreciar en el grafico de 
la figura 8. 




Jo 



M 



Figura 8 Cortante cntico de deslaminacidn. 



Los coeficientes de minoracion, en este caso, son mucho mayores que los empleados para el 
dimensionamiento a flexion, debido a la falta de deformacion plastica de la lamina de material 
compuesto, la variabilidad del sistema de ejecucion, el limitado numero de ensayos disponibles,etc. 
Como ejemplo, la normativa suiza SLA emplea un factor de minoracion de 1.2 en el caso de 
dimensionamiento a flexion y de 1.7 para el dimensionamiento por deslaminacion. 
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4.3 Dimensionamiento del anclaje 

En estructuras reforzadas mediante laminas adheridas, se debe prestar especial atencion a las 
tensiones rasantes en la interfaz hormigon - material compuesto y en los extremos de las laminas. 

Las caractensticas adherentes de estos refuerzos dependen de numerosos factores: tamano y tipo 
del refuerzo, condiciones del sustrato, coeficiente de Poisson, resistencia y confinamiento del 
hormigon, tipo de carga, efectos dependientes del tiempo, cantidad de recubnrniento, preparacion 
del sustrato. tipo y volumen de la fibra y de la matriz, etc. Debido a esta enorme variabilidad, no 
existen critenos de diseno generales, por lo que se precisa un estudio particular de cada refuerzo. 
contrastando con los resultados experimentales aplicables a cada caso. 

En diversos ejemplos presentados en el apartado 6 se realizaron analisis por elementos finitos en 
regimen de tension plana en las vigas reforzadas para estimar la longitud minima de transference de 
rasantes desde las laminas de materials compuesto s hasta la superficie del hormigon, con objeto de 
no superar las tensiones rasantes admisibles (la distribucion tipo de esfuerzos cortante en la interfaz 
hormigon - laminado se detalla en la figura 9, fallando el anclaje cuando la tension rasante, x, 
transmitida al sustrato supera el valor admisible, x cadm ). El modelo constitutive empleado en el 
material compuesto fue elastico lineal y en el adhesivo un modelo de dario isotropico escalar para 
poder representar el comportamiento de union ductil existente entre el laminado compuesto v el 
hormigon( 1 4), ( 1 5), ( 1 6) y ( 1 7). 



T 





Figura 9 Distribucion dc esfuerzos rasantes en la interfaz 
hormigdn - laminado a lo largo de la viga. 



En los laboratonos EMPA se ha investigado y ensayado el anclaje de laminas CFRP. Si bien las 
reg as de diseno practicas no estan disponibles en el momento presente, se recomienda la aplicacion 
e las normas alemanas aprobadas para el calculo de laminas metalicas como criterio conservador 
(basadas en investigaciones de Ramich). En dicha formulacion, adaptada a las laminas CFRP, se 
tienen en cuenta las propiedades mecanicas del laminado, el modulo de corte del adhesivo (factores, 
ambos, importantes en la capacidad de transference de cargas desde el laminado al sustrato de 
hormigon), etc. A contmuacion se expone la expresion de la longitud de anclaje mencionada ( 1 8): 
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I _ {y R z L ,d y 

b ~b L 2 Kt e T k 



> 500 mm 



cony R Z Ld <i LJ A L 



donde: 

l b : 

Yr ; 

b L : 

k: 

Z L,d* 

t,*. 

Xr: 



longitud de anclaje 
coeficiente de seguridad ( 1 .2) 
ancho de la lamina 
tension de referenda 

fuerza a anclar en la seccion de maximo momento 
ancho de una lamina de acero equivalente 
tension adherente ultima, segun ensayos 



4.4 Dimensionamiento a cortante 



El diseno de refuerzos a cortante se realizara por analogia al del calculo de hormigon armado. 

Ante una insuficiencia de armadura de esfuerzo cortante en la estructura existente, incrementos 
de sobrecargas, cambios de uso, etc. se pueden emplear materiales compuestos tanto tipo roving 
como tejido (figura 10); garantizando en todo momento su anclaje y verificando que las tensiones 
rasantes transmitidas al hormigon sean inferiores a los valores admisibles. 

Extsten diversas formulaciones concretas para cada tipo de material ensayado, pudiendo aplicar 
las recomendaciones francesas para los tejidos de carbono con resina epoxfdica, las suizas para 
tejidos similares, tanto de carbono como de vidrio, las japonesas y canadienses (19). 




Figura 10 Aplicacibn de materiales compuestos, tipo 
roving, como refuerzo a cortante en puenles, 
Colombia. 



4.5 Verificacion de estados 1 unites de servicio 

Los Estados Lfmites de Servicio de elementos estructurales reforzados mediante materiales 
compuestos se describen en terminos de limitaciones de flecha y de anchura de fisura. 

Existen numerosas formulaciones empiricas presentes en la literatura para determinar el ancho 
de fisura maximo y la rigidez de este tipo de estructuras, dependiendo de numerosos parametros: 
modulo de deformacion del material compuesto, su seccion y espesor, tipo y espesor de resina 
empleada en la interfaz, etc. 
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Un criterio sencillo y conservador para verificar los Estados Lfmites de Servicio seria la 
aplicacion de los criterios usualmente aceptados para el hormigon armado (los mecanismos 
adherentes que se desarrollan entre el acero interior y el hormigon son menos eficaces que los 
experimentados en los refuerzos mediante laminados compuestos, observandose en los ultimos una 
fisuracion mas distribuida y de menor ancho). 



5 EJECUCION DE REFUERZOS CON MATERIALES COMPUESTOS 



5.1 Preparation del sustrato de hormigon 

El sustrato que recibira las laminas de materiales compuestos debe responder a determinados 
criterios de reception geometricas, mecanicas y fisico - quimicas. 

Con objeto de soportar los esfuerzos rasantes transmitidos al sustrato a traves de la interfaz 
hormigon - laminado compuesto, se recomiendan sustratos de hormigon con resistencia rasante 
admisible de 1.5 MPa (o en su defecto, hormigones con f ck > 20 MPa). 

Los sustratos deben estar limpios de todo producto que pueda afectar a la buena adherencia del 
material compuesto: pintura, grasa, agentes desencofrantes, cuerpos de tipo organico o vegetal. Para 
el caso concreto de los tejidos de fibra de carbono (tipo TEC, etc.), su escasa sensibilidad a la 
humedad permite que los sustratos se presenten ligeramente humedos. 

5.2 Ejecucion de laminados tipo roving 

Preparada la superficie del sustrato, se procede a aplicar un adhesivo adecuado y compatible con 
las laminas empleadas sobre la zona de hormigon que las recibira, con un espesor de 1 mm, 
mediante espatulas. Para la colocation del mismo se deben respetar todas las indicaciones 
especfficas que el fabricante del adhesivo exponga. 

Seguidamente, se colocan las laminas sobre una mesa o tabla y se limpian de todo tipo de 
particulas y suciedades, aplicando el mismo adhesivo anterior sobre las laminas mediante espatula 
con un espesor de 1 - 2 mm. 

Dentro del tiempo abierto del adhesivo, segun la temperatura, se colocan las laminas sobre la 
superficie de hormigon. Utilizando un rodillo, se presiona la lamina contra el adhesivo hasta que 
este rebose por ambos lados de la laming, eliminando las partes sobrantes. 

Por criterios esteticos (las laminas compuestas con fibras de carbono se presentan en color 
negro) se pueden pintar las laminas del color adecuado para hacer menos llamativa su presencia. 

La ejecucion de este tipo de refuerzos se puede apreciar en las figuras 1 1 a 14. 

Es importante resaltar que el plazo de tiempo empleado en todo este proceso descrito no supera 
las dos-tres horas, por lo que lo convierte en un sistema de refuerzo conveniente donde se produzcan 
interrupciones de trafico. 




Figura 11 Viga de horde impactada, antes de reparar. Puente del Drag<5, Barcelona. 
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Figura 12 Impregnation con resina de las ldminas compuestas antes de su 
colocacidn, puente n° 6 del Enlace del Escorial, Madrid. 




Figura 13 Colocacidn de las l«minas sobre la estructura, puente del Drag6, 
Barcelona. 




Figura 14 Elimination del excedente de resina, puente n° 5 del Enlace del Escorial, 
Madrid. 



Los materiales compuestos en el refuerzo de puentes 



89 



5.3 Ejecucion de laminados en forma de tejidos (tipo TFC) 



La ejecucion de este tipo de refuerzos consta de la aplicacion de una primera capa de resina 
epoxidica con medios que permitan la penetracion de la misma a traves de las irregularidades del 
sustrato, asegurando asi una buena impregnacion. 

El tejido de libras de carbono debe aplicarse sobre la capa de resina humedecida, extendiendose 
de un extremo al otro, verificando que no existan pliegues ni que aparezca excesivamente estirada. 
Posteriormente se presiona sobre la resina, permitiendo que la resina impregne el tejido y la 
eliminacion de eventuales burbujas de aire, mediante un rodillo. 

Despues de colocar el tejido se extiende una nueva capa de resina (llamada capa de cierre), 
aportando material epoxidico que permita impregnar el tejido. Esta capa debe aplicarse 
inmediantamente despues de la colocacion del tejido (en ningun caso, despues de la 
pohmeralizacion de la resina). No es determinante que el espesor de la capa de cierre sea 
estrictamente uniforme, pues no afecta a la resistencia del refuerzo. 



6 REALIZACIONES 



A continuacion se describen diversas estructuras reforzadas con laminados de materiales 
compuestos con fibra de carbono y resina epoxidica, tanto en Espana como en Latinoamerica. 

6.1 Puente del Drago, Barcelona 



El puente del Drago, situado en la Avenida de Meridiana sobre el Paseo de Valldaura, 
Barcelona, (una de las vfas de acceso mas importantes a Barcelona desde el norte) recibio el 
impacto de un vehiculo (diciembre, 1.995) en su viga de borde del tablero lado montana - 
Tarragona, provocado por un galibo insuficiente. 

El puente esta constituido por dos tableros de vigas de canto constante. El esquema estatico 
longitudinal de las estructuras se corresponde con el de un tablero continuo biapoyado, con dos 
tramos de luz aproximadamente igual (unos 12 m). Las 15 vigas que conforman el tablero son de 
hormigon armado y seccion transversal rectangular, conectadas por un forjado de compresion de 20 
cm de espesor y mediante traviesas a tercio de la luz. 

Como consecuencia del accidente, la totalidad de la armadura resistente a flexion longitudinal (6 
barras de acero corrugado de 4> 25 mmrde la seccion centro de vano resulto seccionada (figura 1 1). 
Por razones de seguridad estructural y de adecuado comportamiento en servicio, el puente debia ser 
reforzado con urgencia. 

Una vez realizado el analisis estructural del tablero accidentado se pudieron determinar los 
esfuerzos de calculo en la viga de borde danada, constatandose que con la nueva distribucion de 
esfuerzos (provocada por la perdida de rigidez de la viga extrema) no se requerfa el refuerzo de 
vigas adyacentes. Se presentaron diversas alternativas de refuerzo al Ayuntamiento de Barcelona, 
Admimstracion responsable de la gestion del puente del Drago. Entre las diversas opciones 
encontramos: refuerzo mediante recrecido transversal de hormigon en la viga de borde, reponiendo 
el acero con barras de acero corrugado o con acero de pretensar; refuerzo mediante chapa de acero 
encolada con resina epoxidica en la zona inferior de la viga; refuerzo mediante laminas de 
materiales compuestos con libra de carbono en la zona inferior de la viga. 

Debido a los mayores plazos de ejecucion y medios auxiliares requeridos en las dos primeras 
opciones, el Ayuntamiento de Barcelona determino la opcion del refuerzo con materiales 
compuestos como la adecuada para- la reparacion de la estructura. 
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El material compuesto empleado fue el constitufdo por fibras de carbono (tipo C, Toray T 700) y 
resina epoxfdica. Las caracterfsticas del mismo son las siguientes (referidas a la direction 
longitudinal de la fibra): 

Resistencia a traccion: 24.000 Kp/cm2. 

Modulo de deformacion: 1.500.000 Kp/cm2. 

Deformacion ultima: 1 .4 %. 

Contenido volumetrico de fibra: 60 %. 

Densidad: 1 .6 Kp/cm3. 

Resistencia a la temperatura: > 500 °C. 

El diagrama momento - curvatura de la section crftica de la viga reforzada mediante fibras de 
carbono se puede apreciar en la figura 4. Se emplearon laminas con un ancho total de 240 mm. y un 
espesor de 1 .2 mm. El tipo de lamina de material compuesto empleado se puede apreciar en la figura 
2. El espesor de la resina adherente fue de 1 mm. 

En la figura 15 se puede apreciar una description del refuerzo proyectado para el puente del 
Drago y en la figura 16 la situation final de dicha estructura, una vez llevado a cabo el refuerzo. 




Figura 15 Reparacion Puente del Drago. Alzado y secciones. 
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Figura 16 Vista del pucnte del Drago, una vez finalizados los trabajos de refuerzo. 



6.2 Puentes n° 5 y 6 del Enlace del Escorial, Madrid 



Las estructuras n° 5 y n° 6 del Enlace del Escorial se encuentran en el tramo “Villalba - Caloco” 
de la Autopista A6, en Madrid, sufriendo ambas impactos de vehfculos, provocados por galibos muy 
estrictos, en 1.997. 

• Estructura n° 5 

Consta de dos tableros independientes, uno por cada sentido de circulacion, bajo el tronco 
de la A6. Cada uno de los tableros presentan tres vanos isostaticos de unos 12.4 m de luz y 
15.7 m de ancho, formados por ocho vigas de hormigon pretensado enlazadas por un 
forjado superior de hormigon armado de 22 cm. 

El impacto del vehiculo se produjo en el tramo central del tablero que soporta los carriles de 
la A6, concretamente en la viga de borde de la margen izquierda. La viga impactada, junto 
con su viga adyacente, presentaban un canto de 0.83 m. 

Los desperfectos producidos fueron: hormigon del ala inferior y del alma de la viga 
deteriorados, perdida de acero de pretensado (aproximadamente un 80 % del pretensado 
inicial) y de armadura pasiva y revirado del alma de la viga (figura 17). 

• Estructura n° 6 " 

La estructura consta de dos puentes tableros, uno para cada sentido de circulacion, 
constituido cada uno de ellos por tres vanos isostaticos de unos 14 m de luz. 

La viga impactada fue la viga de borde del tramo central, formado por 12 vigas de 
hormigon pretensado, de 0.8 m de canto, enlazadas por un forjado superior de hormigon 
armado de 22 cm. 

Los danos causados fueron similares a los de la estructura n° 5: perdida de acero de 
pretensado y armadura pasiva, hormigon totalmente deteriorado en ala inferior y fisuracion 
en zona superior del alma (figura 1 8). 



Debido al trafico intenso que circunda la zona, se adopto la solucion de refuerzo con materiales 
compuestos, gracias a la rapidez de la ejecucion del mismo. 

Las reparaciones, restituciones y refuerzos propuestos se exponen a continuacion: restitucion del 
hormigon en todas aquellas zonas deterioradas por los impactos, sellado de todas las fisuras 
producidas, recrecido lateral de hormigon en ambos lados de las vigas de borde impactadas con 
objeto de garantizar un comportamiento monolitico de la viga e incremento de la capacidad portante 
de las vigas mediante un refuerzo constituido con laminas de material compuesto con fibra de 
carbono. El tipo de material empleado fue el mismo que el utilizado en el puente del Drago. 
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Figura 17 Deterioros en viga de horde impactada, 
pucnte n° 5 del Enlace del Escorial, 
Madrid. 




Figura 18 Deterioros en viga de horde impactada. 

pucntc n° 6 del Enlace del Escorial. 
Madrid. 



6.3 Realizaciones en Latinoamerica 

• Puente El Consuelo, Colombia 

Puente de hormigon armado de una sola luz de 28 m. Consta de un tablero de 4.6 m de 
ancho y de dos vigas de altura 2.04 m (incluida la losa), con un ancho de alma de 0.32 m. 
Debido a las exploraciones petroliferas que se estan realizando por la zona, se requena 
movilizar sobre la estructura un taladro de perforation de 40 toneladas, por lo que se 
reforzaron ambas vigas a flexion con laminas de fibra de carbono, tipo roving, de 10 cm de 
ancho y 1.2 mm de espesor. La estructura se puede apreciar en la figura 19, donde se 
observa la ligereza de los medios auxiliares empleados para la ejecucion del refuerzo. 

• Puente El Can, Colombia 

Puente de dos tramos de luces desiguales con 7 vigas simplemente apoyadas en pila central, 
la cual consta de 8 columnas unidas en su parte superior por una viga cabezal. La via se 
encuentra local izada en zona de explotacion minera con trafico pesado de camiones. Fue 
necesario reforzar los voladizos de las vigas cabezales en tres puentes de las mismas 
caracterfsticas siaiados en la jnisma via (figura 20), mediante laminas de material 
compuesto de fibra de carbono, tipo roving, de 5 cm de ancho y 1.2 mm de espesor. 




Figura 19 Puente El Consuelo, Colombia. 



Figura 20 Puente El Can, Colombia. 
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• Puente Cocoma, Colombia 

Puente de planta curva constituido por tramos rectos isostaticos. El tramo central esta 
form ado por vigas de hormigon pretensado apoyadas en dos voladizos a traves de mensulas. 
Los voladizos y los tramos compesatorios presentan una seccion viga cajon, mostrando una 
fisuracion inclinada de aproximadamente 45°, prolongandose hasta la losa del tablero. Se 
realizo un refuerzo a cortante mediante laminas de material compuesto de fibra de carbono. 
tip° roving, de 2.5 m de longitud, separadas 90 cm y orientadas perpendicularmente a la 
direccion de las fisuras (figura 10). 



7 CONCLUSIONES 

Con todo lo expresado anteriormente sobre la aplicacion de materiales compuestos en el refuerzo 
de estructuras, podemos concluir con los siguientes puntos: 

• El rango de aplicacion del sistema de refuerzo mediante materiales compuestos es muy 
amplio, pudiendose Ilevar a cabo no solo donde los requerimientos de altas prestaciones y/o 
la maniobrabilidad y rapidez de la ejecucion conslituyan un punto importante, sino, tambien, 
donde el planteamiento sea estrictamente economico, gracias al ahorro que se obtiene por la 
utilizacion de medios auxiliares ligeros con plazos de ejecucion mmirnos (del orden de hasta 
un 25% frente a otras alternativas de refuerzo mas convencionales). La reduccion importante 
de los plazos de ejecucion le convierte en un sistema de refuerzo muy adecuado donde se 
produzcan interrupciones de trafico, etc. 

• A pesar del comportamiento constitutive completamente elastico lineal de los materiales 
compuestos, puede obtenerse una respuesta seccional de la estructura (hormigon + refuerzo) 
ductil si se disena adecuadamente el refuerzo, permitiendo la plastificacion del acero de 
armar antes de que el refuerzo alcance su deformacion ultima. 

• En la fiabilidad de una estructura reforzada con materiales compuestos influye sobremanera 
el tipo de material compuesto, la solicitacion de trabajo, condiciones ambientales, etc., por lo 
que resulta compleja la extrapolacion de los criterios de diseno y de los diferentes 
coeficientes de seguridad para los diferentes materiales compuestos, requiriendose un estudio 
particular para cada material y aplicacion del mismo. 
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